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Abstract 

The relation between limit cycles of planar differential systems and the 
inverse integrating factor was first shown in an article of Giacomini, Llibre 
and Viano appeared in 1996. From that moment on, many research articles 
are devoted to the study of the properties of the inverse integrating factor 
and its relation with limit cycles and their bifurcations. This paper is a 
summary of all the results about this topic. We include a list of references 
together with the corresponding related results aiming at being as much 
exhaustive as possible. The paper is, nonetheless, self-contained in such a 
way that all the main results on the inverse integrating factor are stated and 
a complete overview of the subject is given. Each section contains a different 
issue to which the inverse integrating factor plays a role: the integrability 
problem, relation with Lie symmetries, the center problem, vanishing set of 
an inverse integrating factor, bifurcation of limit cycles from either a period 
annulus or from a monodromic w-limit set and some generalizations. 
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1 The Euler integrating factor 

The method of integrating factors is, in principle, a means for solving ordinary 
differential equations of first order and it is theoretically important. The use of 
integrating factors goes back to Leonhard Euler. 

Let us consider a first order differential equation and write the equation in the 
Pfaffian form 

uj = P{x, y) dy — Q{x, y) dx = . (1) 

We assume that the functions P and Q are of class in a region W C M^. If there 
is a solution of ([1]) which may be expressed in the form H{x, y) = h with H having 
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continuous partial derivatives in lA and with h an arbitrary constant, then it is not 
difficult to see that such an H satisfies the linear partial differential equation 




(2) 



Conversely, every non-constant solution H oi ^ gives also a solution H{x, y) = h 
of ([1]). Thus, solving ([1]) and solving ([2]) are equivalent tasks. 

It is straightforward to show that if Hq{x, y) is a non-constant solution of equa- 
tion ([2]), then all solutions of this equation are of the form F{Hq{x, y)) where F is 
a freely chosen function with continuous derivative. The connection between equa- 
tions ([T]) and ([2]) may be presented also in another form. Suppose that H{x, y) = h 
is any solution of ([1]). Then ([2]) implies 



If we denote the common value of these two ratios by y), then we have 
dH/dy = iiP and dH/dx = —fiQ. This gives to the differential of the function 
H the expression dH{x, y) = y){P{x, y) dy — Q{x, y) dx). Hence, y) is 
called the integrating factor of the given differential equation ([1]) because the left 
hand side of ([T]) turns, when multiplied by y), to be an exact differential. 

Conversely, any integrating factor /i of ([I]), i.e. such that y){P{x, y) dy — 
Q{x, y) dx) is the differential of some function H, is easily seen to determine the 
solutions of the form H{x, y) = h of iQ. Altogether, solving the differential 
equation ([1]) is equivalent to finding an integrating factor of the equation. 

When an integrating factor /i of ([T]) is available, the function H can be obtained 
from the line integral 



along any curve connecting an arbitrarily chosen point (xq, yo) and the point (x, y) 
in the region U. We remark that this line integral might not be well-defined if the 
region U is not simply-connected. When we know an integrating factor yU of ([T]), 
we have a first integral well-defined in each simply-connected subcomponent of the 
region U. 

2 The inverse integrating factor 

Let us consider a real planar autonomous differential system 



dH/dy 
P 



dH/dx 




X = P{x,y) 



y = Q{x,y) 



(3) 
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where P{x,y) and Q{x,y) are of class C^iU) and W C is an open set. The dot 
denotes derivation with respect to the independent variable t usually called time, 
that is ' = 4. 

at 

As usual, we associate to system ([3]) the vector field X = P{x, y)dx + Q{x, y)dy. 
Notice that the ordinary differential equation a; = given in ^ is just the differ- 
ential equation of the orbits of system ([3]). 

Definition 1 A function V : W — M said to he an inverse integrating factor of 
system ([3]) if it is of class C^(U), it is not locally null and it satisfies the following 
partial differential equation: 

.dV{x,y) dV{x,y) fdP{x,y) dQ{x,y)\ 

In short notation, an inverse integrating factor V of system satisfies XV = 
VdivX, where divA* = §f + ^ stands for the divergence of the vector field X. 

Of course, the computation of an inverse integrating factor for a concrete system 
is a delicate matter whose difficulty is comparable to solving the system itself. 

If V is an inverse integrating factor of a vector field X, then the zero set of 
V, V^~^(0) := {{x,y) \ V{x,y) = 0}, is composed of trajectories of X. For by the 
equation (jlj) that defines V, X is orthogonal to the gradient vector field W along 
the zero set of V. 

The name "inverse integrating factor" arises from the fact that if V solves 
equation (jll), then its reciprocal 1/V is an integrating factor for X onU \ ^^^(0). 

3 Local nontrivial Lie symmetries and inverse in- 
tegrating factors 

Roughly speaking, a symmetry group of a system of differential equations is a con- 
tinuous group which transforms solutions of the system to other solutions. Simple 
typical examples are groups of translations, rotations and scalings, but these cer- 
tainly do not exhaust the range of possibilities. Once one has determined the 
symmetry group of a system of differential equations, a number of applications 
become available. 

More precisely, a symmetry of system Qin U, where W C is an open set, is a 
1-parameter Lie group of diffeomorphisms $e acting in U that maps the set of orbits 
of ([3]) into itself. When y) = {x{x, y; e),y{x, y; e)), the symmetry condition of 
([2D reads for x = P{x,y), y = Q{x,y) for all e close to zero. Let the C^{U) vector 
field y = ^{x, y)dx + rj{x, y)dy be the infinitesimal generator of the 1-parameter Lie 
group that is, x{x, y;e) = x + e^(x, y) + O(e^), y{x, y;e) = y + eri{x, y) + O(e^). 
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Denoting by X = P{x,y)dx + Q{x,y)dy the vector field associated to system ([3]), 
it is well known that a characterization of the Lie symmetries of ([3]) is given by 
the relation [Af, 3^] = fi{x, y)X for certain scalar function /i : W — > M. In this 
expression we have used the Lie bracket of two C^-vector fields X and y defined as 
[A", 3^] := Xy — yX. Using coordinates we have 

r;f 3;] = (p^ _ + _ r^^] + (p^ _ + _ r^^] Q (5) 
\ dx dx dy dy J ^ \ dx dx dy dy J ^ 

When beginning students first encounter ordinary differential equations, they are 
presented with a variety of special techniques designed to solve certain particular 
types of equations, such as separable, homogeneous or exact. Indeed, this was the 
state of the art around the middle of the nineteenth century, when Sofus Lie made 
the profound discovery that these special methods were, in fact, all special cases of 
a general integration procedure based on the invariance of the differential equation 
under a continuous group of symmetries. This observation at once unified and 
significantly extended the available integration techniques. 



Differential Equation 


Lie Symmetry 


dy/dx = f{x)g{y) 


y = 9iy)dy 


dy/dx = f{ax + by) 


y = bd^ + ady 


dy/dx - 


y = ydx- xdy 


dy/dx = f{y/x) 


y = xdx + ydy 


dy/dx = P{x)y + Q{x) 


y = exp (J P{x)dx) dy 


dy/dx = P{x)y + Q{x)y'^ 


y = y"' exp {1 — n) J P{x)dx dy 



Consider now a vector field X = P{x,y)dx + Q{x,y)dy defined in an open 
connected subset C M^. In the case of a single first order ordinary differential 
equation dy/dx = Q{x, y)/P{x, y), the Lie symmetries method provides by quadra- 
ture an explicit formula for the general solution. In fact, one can easily see that if 
we know a Lie symmetry in U with infinitesimal generator y = C,{x, y)dx + r]{x, y)dy 
then we construct an inverse integrating factor V = det{X, y} = Prj — QC, defined 
in U, but the converse is not always true. To see that, assume now the existence 
of an inverse integrating factor of A" in a simply connected domain U and we 
look for an infinitesimal generator y = C,{x,y)dx + ri{x,y)dy of a Lie symmetry of 
X well defined in U. We recall that a singular point p ^ U of A" is called weak if 
divX{p) = 0. If there is no weak singularity of A" in W, then we can do at least 
one of the following constructions: 

(i) Prescribe the function C,{x,y) and solve r]{x,y) from V = Prj — Q^. 
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(ii) Prescribe the function r]{x,y) and solve C,{x,y) from V = Prj — Q^. 



(iii) Take the rescaled hamiltonian vector field 



1 dV dV 
dwX dy ox 



defined in U\{{x, y) eU : AivX = 0}. 

Therefore, the equivalence between inverse integrating factors and Lie symmetries 
for planar vector fields X is not true, in general, in neighborhoods of weak singular 
points of X. Of course, some special situations can appear giving the equivalence 
when X possesses an analytic first integral in these neighborhoods as the nonde- 
generate center singular point shows. 



Importance of inverse integrating factors arises from the fact that the differential 
1-form uj/V = {P dy~Q dx)/V is closed {d{uj /V) = 0) in U\V-^{0). Then in the 
case in which U\V~^{0) is simply-connected, the 1-form u/V is exact {w/V = dH), 
and therefore a first integral H{x,y) of the differential equation is immediately 
constructed. As a consequence, the vector field X = P{x, y)dx + Q{x, y)dy is topo- 
logically equivalent, in U, to the hamiltonian vector field X /V = — 

Making a pause in this exposition we now present an example. Let us consider 
the following cubic system 

X = P{x,y) = -y - x{x^ + y'^ -I) , y = Q{x,y) = x - y{x^ + y"^ - I) . (6) 

An inverse integrating factor for system (jS]) is given by V{x,y) = (x^ + + 
y"^ — 1). Associated to him one has the first integral 

Hix, y) = exp {2 arctan (^) } , 

which is not continuous in (0,0). On the other hand, since the polar form of the 
system is r = 2r^(r^ — 1), = 1 it is easy to check that the unit circle x^+y^ — 1 = 
is the unique limit cycle of system (EI). Let X = P{x, y)d/dx + Q{x, y)d/dy be the 
vector field associated with system ([6]). From the symmetries point of view, since 
3^ = ydx — xdy satisfies [.Y, 3^] = we have that y is the infinitesimal generator 
of a Lie group admitted by system ([6]) which is just the 5*0(2) rotation group 
X = X cos e — ?/ sin e, y = xsme + y cos e. Hence V{x, y) = det{X, y} is an inverse 
integrating factor of system (Q. Notice that the only common integral curves for 
the vector fields X and y are included in V~^{0) and are just the separatrices of 
X. This behavior will be explained in future sections. 
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4 On the integrability problem 



The integrability problem is mainly related to planar polynomial differential sys- 
tems of the form 

x = P{x,y), y = Q{x,y), (7) 

where P{x,y),Q{x,y) G are coprime polynomials, that is, there is no non- 

constant polynomial which divides both P and Q. We call d the maximum degree 
of P and Q and we say that system ([7]) is of degree d. When d = 2, we say that 
(|7j) is a quadratic system. 

If p is a point such that P{p) = Q{p) = 0, then we say that p is a singular point 
of system ([7]). 

As we have already defined in Section [1], a function H -.U —>-M. such that it 
is constant on each trajectory of ([7]) and it is not locally constant is called a first 
integral of system ([7]) of class j defined on W C M?. The equation H{x, y) = h for a 
fixed h gives a set of trajectories of the system, but in an implicit way. When 
j > 1, these conditions are equivalent to P(x,y)^ + Q{x,y)^ = and H not 
locally constant. The problem of finding such a first integral and the functional 
class it must belong to is what we call the integrability problem. 

To find an integrating factor or an inverse integrating factor for system ([7]) is 
closely related to finding a first integral for it. When considering the integrability 
problem we are also addressed to study whether an (inverse) integrating factor 
belongs to a certain given class of functions. 

When a first integral H of system ([7]) is known, all the orbits of the system are 
contained in its domain of definition are given by the level sets H{x, y) = h. Thus, 
a natural strategy is to look for the determination of some of the orbits of the 
system and try to build a first integral with them. In particular, and since system 
([7j) is polynomial, those orbits which are algebraic will be of special interest. 

An invariant curve is a curve given by f{x,y) =0, where / : W C — >• M is a 
function in the open set U, non locally constant and such that there exists a 
function in U, denoted by k{x,y) and called cofactor, which satisfies: 

df df 
P{x,y)—{x,y) + Q{x,y) —{x,y) = k{x,y) f{x,y), (8) 

for all {x,y) G U. The notion of invariant curve was first introduced in ^53j. The 
identity (IHl) can be rewritten by Xf = kf. We recall that Xf denotes the scalar 
product of the vector field X and the gradient vector V/ related to f{x,y), that 
is, Vf{x,y) = {^{x,y), ^{x,y)). We will denote by f or by / the function Xf 
once evaluated on a solution of system ([7]). In case f{x,y) = defines a curve in 
the real plane, this definition implies that the function Xf is equal to zero on the 
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points such that f{x,y) = 0. In the article [53] an invariant curve is defined as 
a function f{x,y) defined in the open set U C M?, such that, the function Xf 
is zero in all the points {{x,y) G U \ f{x,y) = 0}. We notice that our definition 
of invariant curve is a particular case of the previous one but, for the sake of 
our results, the cofactor is very important and that's why we always assume its 
existence. 

When the cofactor k{x, y) is a polynomial, we say that /(x, y) = is an invariant 
curve with polynomial cofactor. We only admit invariant curves with polynomial 
cofactor of degree lower or equal than d — 1, that is degk{x,y) < d — 1, where d 
is the degree of system ([7]). 

The notion of invariant curve is a generalization of the notion of invariant 
algebraic curve. An invariant algebraic curve is an algebraic curve f{x,y) = 0, 
where f{x, y) G C[x, y], which is invariant by the fiow of system ([7]). This condition 
equals to Xf = kf, where the cofactor of an invariant algebraic curve is always a 
polynomial of degree deg k{x, y) < d — 1. 

We cite [78l[99|IT00] as compendiums of the results on invariant algebraic curves. 
For instance, in [7S], it is shown that if f{x, y) = and g{x, y) = are two invariant 
algebraic curves of system ([7]) with cofactors kf{x,y) and kg{x,y), respectively, 
then the product of the two polynomials gives rise to the curve {fg){x,y) = 
which is also an invariant algebraic curve of system ([7]) and whose cofactor is 
kfix,y) + kg{x,y). 

In order to state the known results of integrability using invariant algebraic 
curves, we need to consider complex algebraic curves f{x,y) = 0, where f{x,y) G 
C[x, y]. Since system ((Tj) is defined by real polynomials, if f{x, = is an invariant 
algebraic curve with cofactor k{x,y), then its conjugate f{x,y) = is also an 
invariant algebraic curve with cofactor k{x,y). Hence, its product f{x,y)f{x,y) G 
M.[x, y] gives rise to a real invariant algebraic curve with a real cofactor k{x, y) + 
k{x,y). For a sake of simplicity, we consider invariant algebraic curves defined by 
polynomials in C[a:, ?/], although we always keep in mind the previous observation. 
In M^, the curve given by f{x,y) = 0, where f{x,y) is a real function, may only 
contain a finite number of isolated singular points or be the null set. 

An algebraic curve f{x, ?/) = is called irreducible when f{x, y) is an irreducible 
polynomial in the ring C[a:, ?/]. We can assume, without loss of generality, that 
f{x,y) is an irreducible polynomial in C[x,y], because if f{x,y) is reducible, then 
all its proper factors give rise to invariant algebraic curves. Given an algebraic curve 
f{x,y) = 0, we can always assume that the polynomial f{x,y) has no multiple 
factors, that is, its decomposition in the ring C[a;, is of the form f{x,y) = 
fi{x, y)f2{x, y) . . . fiix, y), where fi{x, y) are irreducible polynomials and fi{x, y) ^ 
cfj{x, y) iii ^ j and for any c G C. The assumption that given an algebraic curve 
f{x,y) = 0, the polynomial f{x,y) has no multiple factors is mainly used to 
ensure that we do not consider "false" singular points. If p is a point such that 
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f{p) = and V fij)) = 0, and f{x,y) has no multiple factors, then p is a singular 
point of the curve f{x,y) = 0. But, if f{x,y) has multiple factors, for instance, 
f{x,y) = fi{x,yY where fi{x,y) is an irreducible polynomial in C[x,y], then all 
the points of the curve {p\ fi{p) = 0} satisfy the property that f{p) = and 
V/(p) = although they are not all singular points. 

We recall that if p is a singular point of an invariant algebraic curve f{x, y) = 
of a system ([7]) , then p is a singular point of the system. Given an algebraic curve 
f{x,y) = 0, we will always assume that the decomposition of f{x,y) in the ring 
C[x, ?/] has no multiple factors. We want to generalize this property to invariant 
curves, that's why we will always assume that, given an invariant curve f{x,y) = 0, 
if p G W is such that f{p) = and V f{p) = 0, then p is a singular point of system 
([7j) . This technical hypothesis generalizes the notion of not having multiple factors 
for algebraic curves. In [23], a set of necessary conditions for a system ([7]) to have 
an irreducible invariant algebraic curve is given. 



Invariant algebraic curves are the main objects used in the Darboux theory 
of integrability. In [38J, G. Darboux gives a method for finding an explicit first 
integral for a system ([7]) in case that d(d + l)/2 + 1 different irreducible invariant 
algebraic curves are known, where d is the degree of the system. In this case, a 
first integral of the form H = /^V2 ^ • • • fs% where each fi{x, y) = is an invariant 
algebraic curve for system ([7]) and Aj G C not all of them null, for i = 1,2, ... ,s, 
s G N, can be constructed. The functions of this type are called Darboux functions. 

As we have already stated, given an invariant algebraic curve f{x, y) = whose 
imaginary part is not null, then its conjugate is also an invariant algebraic curve. 
Moreover, as system ([7]) is real, if f{x, y) appears in the expression of a first integral 
of the form given by Darboux with exponent A, then f{x, y) appears in the same 
expression with exponent A. We call Re/ the real part of the polynomial / and by 
Im/ its imaginary part. Analogously, let us call ReA the real part of the complex 
number A and by ImA its imaginary part. We call i = V— 1 and we use the following 
formula for complex numbers: 



arctan(2;) = log 



1 



iz 



iz 



i/2 



zeC, 



to show that 



(Re/ + Im/ 1)^^=^+!-^' (Re/ - Im/ if-^-''^^' 
{(Reff + (Imff)^"^ exp |-2ImA arctan 



Im/ 
Re/ 



We deduce that the product f{x,y)^f{x,y)^ is a real function and so it is any 
Darboux function H = f^'f^^ . . . 



8 



We have that the Darboux function H can be defined in the open set \ S, 
where S = {(x, ?/) G | (/i ■ /2 ■ • • • ■ fr){x, y) = 0}. We remark that, particularly, 
if Aj G Z , Vi = 1,2, ... ,r, if is a rational first integral for system ([7j). In this 
sense J. P. Jouanoulou [76], showed that if at least d(d + 1) + 2 different irreducible 
invariant algebraic curves are known, then there exists a rational first integral. 

The main fact used to prove Darboux's theorem (and Jouanoulou's improve- 
ment) is that the cofactor corresponding to each invariant algebraic curve is a 
polynomial of degree < d — 1. Invariant curves with polynomial cofactor can also 
be used in order to find a first integral for the system. This observation enables a 
generalization of the Darboux's theory which is given in [52], where, for instance, 
non-algebraic invariant curves with an algebraic cofactor for a polynomial system 
of degree 4 are presented. In [50], other examples are given of such invariant 
curves with polynomial cofactor for some families of systems and the way they 
are used to construct explicit first integrals and inverse integrating factors for the 
corresponding systems. As a continuation of [60j, in [61] we study when a planar 
differential system polynomial in one variable linearizes in the sense that it has an 
inverse integrating factor which can be constructed by means of the solutions of 
linear differential equations and we describe some families of differential systems 
which are Darboux integrable and whose inverse integrating factor is constructed 
using the solutions of a second-order linear differential equation defining a family 
of orthogonal polynomials. 

Some generalizations of the classical Darboux theory of integrability may be 
found in the literature. For instance, independent singular points can be taken into 
account to reduce the number of invariant algebraic curves necessary to ensure 
the Darboux integrability of the system, see |29j. A good summary of many of 
these generalizations can be found in [90] and a survey on the integrability of two- 
dimensional systems can be found in ^2T\. One of the most important definitions 
in this sense is the notion of exponential factor which is given by C. Christopher in 
[31| . when he studies the multiplicity of an invariant algebraic curve. The notion 
of exponential factor is a particular case of invariant curve for system ([7j). Given 
two coprime polynomials h,g E y], the function e'^/^ is called an exponential 
factor for system ([7]) if for some polynomial k of degree at most d — 1 , where d is 
the degree of the system, the following relation is fulfilled: 



The next proposition, proved in [31], gives the relationship between the notion 
of invariant algebraic curve and exponential factor. 




As before, we say that k{x, y) is the cofactor of the exponential factor e'^/^. 



9 



Proposition 2 [31| If F = e^^^ is an exponential factor and g is not a constant, 

dx^^ dy 



then g = is an invariant algebraic curve, and h satisfies the equation P^^^+Q^ 



hkg + g kp where kg and kp are the cofactors of g and F, respectively. 

The notion of exponential factor is very important in the Darboux theory of in- 
tegrability since it does not only allow the construction of first integrals following 
the same method described by Darboux, but it also explains the meaning of the 
multiplicity of an invariant algebraic curve in relation with the differential system 
(I7j). A complete work on this subject can be found in [34|. 

In the same way as with invariant algebraic curves, given an exponential factor 
F = exp{h/g}, since system ([7]) is a real system, there is no lack of generality in 
considering that h{x, y),g{x, y) G M[x, y]. li F = exp{h/g} is an exponential factor 
with non-null imaginary part, then its complex conjugate, F = exp{h/g} is also an 
exponential factor, as it can be easily checked by its defining equation. Moreover, 
the product F F = exp{h/g + h/g} is a real exponential factor with a real cofactor. 

Since the notion of exponential factor is the most current generalization in the 
Darboux theory of integrability, any function of the form: 

ft ft ■ ■ ■ fr'- f-p V" f-p y ' ■ • • fexp f y* , (9) 



\ \9i'JJ V \92 J J V \9e / / 

where r, £ G N, fi{x, y) = {1 < i < r) and gj{x, y) = {1 < j < i) are invariant 
algebraic curves of system ([7]), hj{x,y) (1 < j < i) are polynomials in C[x,y], Aj 
{1 < i < r) and /ij (1 < j < i) are complex numbers and fij {1 < j < i) are 
non-negative integers, is called a (generalized) Darboux function. 

Let us present a short survey about the Darboux method and its improvements. 
Let us recall that a singular point (xq, yo) of system ([7]) is called weak if the diver- 
gence, divX, of system ([7j) at (xo,?/o) is zero. We recall that X denotes the vector 
field associated to system ([7j). We denote by Cd-i[x,y] the set of polynomials in 
C[x, y] of degree lower than d. We say that s points {xk, yu) G C^, /c = 1, 2, . . . , s, 
are independent with respect to Cd-i[a;,i/] if the intersection of the s hyperplanes 

*+J=0 ) fc=l,2,, 

is a linear subspace of C'^('^+^)/^ of dimension d(d + l)/2 — s > 0. 

The main results about the Darboux method and its improvements are sum- 
marized in the following theorem, which can be found in |80], see also |90j . 
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Theorem 3 Suppose that a polynomial differential system ([7j) of degree d admits 
r irreducible invariant algebraic curves = with cofactors Ki for i = 1,2, . . . ,r; 
i exponential factors exp{hj/ gj^) with cofactors Lj for j = 1,2, . . . , i; and s inde- 
pendent singular points {xk, yk) such that fi{xk, yu) 7^ for i = 1,2, ... ,r and for 
k = 1,2, ... ,s. Moreover, the irreducible factors of the polynomials gj are some 

f^'s. 

(a) There exist Xi, fij G C not all zero such that ^[^j^ + Y^^j=if^j^j — ^> 
and only if the (multivalued) function <^ is a first integral of system ([7]). 

(b) If r + i + s = [d(d + l)/2] + 1, then there exist Xi, fij G C not all zero such 
that X;i=i XiKi + Y!'j=i l^jLj = 0. 

(c) Ifr + i + s > [d(d + l)/2] +2, then system has a rational first integral, and 
consequently all trajectories of the system are contained in invariant algebraic 
curves. 

(d) There exist Xi, fij G C not all zero such that J2l=i ^t^i + X]j=i f^j^j = ^iv X 
if and only if the function ([9]) is an inverse integrating factor of system (JTj). 

(e) If r + i + s = d(d + l)/2 and s independent singular points are weak, then 
the function ([9]) for convenient Xi, fij G C not all zero is a first integral or an 
inverse integrating factor of system (171). 

Introducing the notion of muhiphcity of the invariant algebraic hypersurfaces 
of a polynomial vector field in C^, the results of Darboux integrability theory of 
Theorem [3] have been generalized to systems in C^, where n>2, see [82] and the 
references therein. 

An improvement of the previous Darboux theorem is presented in [18] when 
the system has a center. As usual [gJ means the integer part of the real number 

Theorem 4 [T8] Consider a polynomial system ([7]) of degree d, with a center at 
the origin and with an arbitrary linear part. Suppose that this system admits d(d + 
l)/2 — L(d + 1)/2J invariant algebraic curves or exponential factors. Then this 
system has a Darboux inverse integrating factor. 

In the following section we present several relations between the existence of an 
inverse integrating factor and the center problem. 

We recall that the integrability problem consists in finding the class of functions 
a first integral of a given system ([7]) must belong to. We have system ([7]) defined 
in a certain class of functions, in this case, the polynomials with real coefficients 
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]R[x, y], and we consider the problem whether there is a first integral in another, 
possibly larger, class. For instance in [22], H. Poincare stated the problem of 
determining when a system ([7]) has a rational first integral. The works of M.J. 
Prelle and M.F. Singer and M.F. Singer [IQl] go on this direction since they 
give a characterization of when a polynomial system ([7]) has an elementary or 
a Liouvillian first integral. An important fact of their results is that invariant 
algebraic curves play a distinguished role in this characterization. Moreover, this 
characterization is expressed in terms of the inverse integrating factor. 

Roughly speaking, an elementary function is a function constructed from ratio- 
nal functions by using algebraic operations, composition and exponentials, applied 
a finite number of times, and a Liouvillian function is a function constructed from 
rational functions by using algebraic operations, composition, exponentials and in- 
tegration, applied a finite number of times. A precise definition of these classes of 
functions is given in [9^ 1101] . We are mainly concerned with Liouvillian functions 
but we will state some results related to integration of a system ([7]) by means of 
elementary functions. 

We recall that C(x, y) denotes the quotient field associated to the ring of poly- 
nomials with complex coefficients, that is, C{x, y) is the field of rational functions 
with complex coefficients. 

Theorem 5 [I94j // the system (|7]) has an elementary first integral, then there 
exist Wo, wi, . . . ,Wn algebraic over the field C(x, y) and ci, C2, . . . , in C such that 
the elementary function 



H = wo + ^Ci\n{wi) (10) 



=1 



is a first integral of system ([7]) . 



The existence of an elementary first integral is intimately related to the existence 
of an algebraic inverse integrating factor, as the following result shows. 

Theorem 6 [I94j // the system ([7]) has an elementary first integral, then there is 



an inverse integrating factor of the form 

^ ( Aix,y) \ 
\B{x,y)J ' 

where A,B & C[x, y] and N is an integer number. 

The paper [22] is devoted to study which is the form of the inverse integrating 
factor of a polynomial planar system ([7]) with a Darboux first integral H of the 
form Q. This work is an improvement of the results of Prelle and Singer in [HI] 
where it is shown that these Darboux integrable vector fields have a rational inverse 
integrating factor (see Theorem 7 of [94J). In [22], another proof of this result is 
presented. 
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Theorem 7 [22| If the system ([7]) has a (generalized) Darboux first integral of 
the form then there is a rational inverse integrating factor, that is, an inverse 
integrating factor of the form: 



B{x,y)' 

where A, B E C[x, y]. 

Unfortunately, not all the elementary functions of the form ffTOl) are of (generalized) 
Darboux type. That's why, we can find systems with an elementary first integral 
and without a rational inverse integrating factor. The following example is of this 
type. The system appears in the works of Jean Moulin-Ollagnier [851 ES], although 
he does not give an explicit expression for the first integral. The Lotka-Volterra 
system: 

i = x(l-^ + y^, y = 2/ (-3 + I -y) , (11) 

has the irreducible invariant algebraic curves x = 0, y = and f{x,y) = 0, where 
f{x,y) := (x — 2)^ — 2xy. Applying the results described in [23], it can be shown 
that this system has no other irreducible invariant algebraic curve. The function 
V{x, y) = x~^^'^y^^'^ f{x, y) is the only algebraic inverse integrating factor of system 
( ITTj) (modulus multiplication by non null constants). Since there is no rational 
inverse integrating factor, we deduce, by Theorem [TJ that there is no (generalized) 
Darboux first integral. An elementary first integral for this system, which is of the 
form (fTOi) . is given by: 

H{x, y) := \p2\fx^ + ln(x — 2 + \/2^/x^) — ln(x — 2 — \p2\fx^) . 

We remark that both Theorems [6] and [7] give a necessary condition to have an 
elementary or (generalized) Darboux, respectively, first integral. The reciprocals 
to the statements of Theorems [6] and [7] are not true. A result to clarify the easiest 
functional class of the first integral once we know the inverse integrating factor 
appears in [16], see also [19|, where the following theorem is stated: 

Theorem 8 [T6| // the system ([7]) has a rational inverse integrating factor, then 
the system has a (generalized) Darboux first integral. 

In any case, the following Theorem [9] ensures that given an algebraic inverse 
integrating factor, there is a Liouvillian first integral. The Liouvillian class of func- 
tions contains the rational, algebraic, Darboux and elementary classes of functions. 

M.F. Singer shows in [101] the characterization of the existence of a Liouvillian 
first integral for a system ([7j) by means of its invariant algebraic curves. 
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Theorem 9 [lUlj System ([7]) has a Liouvillian first integral if, and only if, there 



is an inverse integrating factor of the form V = exp | fl^xo^y^^V^ > where r] is a 
rational 1-form such that drj = 0. 

We recall that when 1-form t] is such that dr] = 0, we say that it is closed and if 
there exists a function such that rj = dip, we say that rj is exact. 

Taking into account Theorem[9], C. Christopher in [32] gives the following result, 
which makes precise the form of the inverse integrating factor. 



Theorem 10 [32| // the system ([7]) has an inverse integrating factor of the form 



exp I /(I^q'I^q)^I ) where rj is a rational 1-form such that drj = 0, then there exists 
an inverse integrating factor of system ([7]) of the form 

V = exp{D/E}l[cl% 

where D, E and the Ci are polynomials in x and y and /j G C. 

We notice that Cj = are invariant algebraic curves and exp{D/E} is an expo- 
nential factor for system ([7j). In fact, since system ([7]) is a real system, we can 
assume, without loss of generality, that is a real function. 

Theorem [TU] states that the search for Liouvillian first integrals can be reduced 
to the search of invariant algebraic curves and exponential factors. Therefore, if 
we characterize the possible cofactors, we have the invariant algebraic curves of a 
system and, hence, its Liouvillian or non Liouvillian integrability. 

Several works study the relation between the existence of invariant algebraic 
curves and the integrability of the system. The existence of an inverse integrat- 
ing factor and the functional class it belongs to is crucial in the resolution of the 
integrability problem, as Theorems [TJ M and [10] show. A number high enough of 
invariant algebraic curves of system ([7]) implies its integrability in one of the ratio- 
nal, elementary or Liouvillian class, due to Darboux's theorem and Jouanoulou's 
improvement, see also Theorems [T6] and [T71 The degree of an invariant algebraic 
curve is not necessarily related with the integrability class of the system, see |2Sl ES] 
and the references therein. 

We conclude this part with a theorem that summarizes some relations between 
inverse integrating factors and first integrals of polynomial vector fields. 

Theorem 11 Let X be a planar polynomial vector field. 

(i) If X has a Liouvillian first integral, then it has a Darboux inverse integrating 
factor. 
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(ii) // X has a Darhoux first integral, then it has a rational inverse integrating 
factor. 



(iii) If X has a polynomial first integral then it has a polynomial inverse integrating 
factor. 

Statement (i) of Theorem [TT] was proved in |1U1] and [32] and statements (ii) and 
(iii) in [22]. 

Another problem related with the inverse integrating factor and the integra- 
bility problem is an inverse problem: given a function V{x,y), the question is to 
find (all the) planar differential systems with V{x, y) as inverse integrating factor. 
In the case of searching for a Darboux inverse integrating factor, a very exhaus- 
tive approach to this problem is given in [331 EOl [90]. The main result of [33] 
establishes, under two generic conditions, all the planar polynomial differential 
systems with an inverse integrating factor of the form V{x.,y) = fi^f2^ ■ ■ ■ fs% 
where each fi{x,y) = is an invariant algebraic curve of the system and A, G C, 
for i = 1, 2, . . . , s, s G N. We do not reproduce the main result of [33j because a 
lot of notation would need to be introduced. 

In [T7] another method to construct systems with a given inverse integrating 
factor is described. In fact, in 1997 the function V was named null divergence 
factor. This method is a generalization of the classical Darboux method to generate 
integrable systems. One of the main results in this paper is the following one. 

Theorem 12 ^ Let X.-, = Pi{x,y)dx + Qi{x,y)dy, with i = l,2,...,n, be 
vector fields defined in an open subset lA C , which have inverse integrating 
factors Vi{x,y) , respectively. Then, the vector field X = P{x,y)dx + Q{x,y)dy with 

p = ^of^ + E^^f n ^.V- 

y i=i \j=i,j^i J 

where Aj are arbitrary real numbers fori = 0, 1, 2, . . . , n, has the inverse integrating 

n 

factor V{x,y) given byV{x^y) = ]^^(a;, y). 

1=1 

Indeed, if two systems have the same inverse integrating factor, a more general 
system which has such inverse integrating factor can be constructed, as it is shown 
in the following proposition. 
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Proposition 13 [T7] Let Xi = Pi{x,y)dx + Qi{x,y)dy with i = 1,2, be two 
vector fields defined in an open subset U C M^, which have the same inverse inte- 
grating factor V{x, y) . Then, the vector field X\ + \X2 has also the function V{x, y) 
as an inverse integrating factor, for arbitrary values of the real parameter A. 

This proposition establislies tliat tlie set of vector fields with the same inverse 
integrating factor forms a M vector space. 

A polynomial inverse integrating factor allows the study of the dynamics of 
system ([7]), because a first integral can be computed, but it is not so involving as 
looking for a polynomial first integral. Indeed, once the degree of a polynomial 
inverse integrating factor is fixed, by an ansatz for instance, the problem of looking 
for it is reduced to a system of linear equations on its coefficients. Many authors 
have used this idea to find families of planar polynomial differential systems of 
the form ([7j) for which all the dynamics can be determined through an inverse 
integrating factor. 

In [20], necessary conditions for a planar polynomial vector field to have a 
polynomial inverse integrating factor are obtained, see also ^\ . All the quadratic 
systems with a polynomial inverse integrating factor are determined in [36j and all 
the quadratic systems with a polynomial first integral are given in [H]. 

In [9] all polynomial first integrals of the non-homogeneous two-dimensional 
Lotka-Volterra system of ordinary differential equations are determined and the 
role of polynomial inverse integrating factors is emphasized. Indeed, new first in- 
tegrals of this class of systems having a polynomial inverse integrating factor is 
presented. The Liouvillian integrability of Lotka-Volterra systems has been stud- 
ied in [861 E]. 

In the work [LSI, planar differential systems of the form ([7]) and defined by 
the sum of homogeneous vector fields are studied. In particular systems with 
degenerate infinity are taken into account. Let us denote by Pd{x,y) and Qd{x,y) 
the terms of the highest degree d in system ([7]). We say that system ([7]) is of 
degenerate infinity if xQd{x, y) —yPd{x, y) = 0. We remark that when a system ([7]) 
with degenerate infinity is embedded into a compact space (either by the Poincare 
compactification into an sphere or when it is embedded in the complex projective 
plane) the line at infinity is filled with singular points. 

We recall that a real function H{x,y) is said to be p-degree homogeneous if 
H{Xx,Xy) = \^ H{x,y) for all (x, in the domain of definition of H{x,y) and 
for all A G M, where p E 

One of the main results in [13] is the following one. 
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Theorem 14 [13] Let us consider the following planar polynomial differential sys- 
tem 

X = Pn{,x,y) + xAA-i{x,y), y = Q„(x,?/) + yAd-i{x,y), (12) 

where Pnix, y) and Qnix, y) are homogeneous real polynomials of degree n, Ad-i(x, y) 
is a real homogeneous polynomial of degree d — 1 and d > n > 1. Let us also con- 
sider the related homogeneous polynomial differential system: 

X = Pnix,y), y = Qn{x,y). (13) 

Then, the following statements hold. 

(a) IfH{x, y) is ap-degree homogeneous first integral of system f|T3l) . then H{x, y) 
is a particular solution of system (fT2l) . 

(b) The homogeneous function Vn+i{x,y) := xQn{x,y) — yPn{x,y) is an inverse 
integrating factor of system (fT3!) . 

(c) The homogeneous function Vn+i{x,y) := xQn{x,y)—yPn{x,y) is a particular 
solution of system ffT2|) . 

(d) If H{x,y) is ap-degree homogeneous first integral of system (fT3|) . Then, the 
function 

d — n 

[xQnix, y) - yPnix, y)) H{x, y)~ 
is a (generalized) Darhoux inverse integrating factor of system ( fT2l) . 

The degree of a polynomial inverse integrating factor in relation with the de- 
gree d of the system can be bounded under certain conditions. The conditions 
established in the following result come from the embedding of a planar vector 
field in CP^, see [23] for the complete definition of this embedding. Consider the 
polynomial differential system (JTj) with P and Q coprime polynomials of maximum 
degree d. Extending system ([71) to a differential equation in the complex projective 
plane CP^, a point (Xq : Iq ^ 0) G CP^ is termed infinite singular point of system 
([7j) if (Xo, Yq) G is a root of the homogeneous polynomial yPd{x, y) —xQd{x, y). 
Here Pa and Qa denote the highest homogeneous components of P and Q of de- 
gree d. Following Seidenberg, a singular point is called simple if the eigenvalues 
A,/i G C associated to its linear part satisfy A 7^ 7^ and A//i ^ Q"*", where Q"*" 
stands for the positive rational numbers. Given a polynomial P(x, y) of degree n, 
we denote by F{X, Y, Z) its projectivization, that is, the homogeneous polynomial 
P(X, y, Z) = Z"'F{X/Z, Y/Z). Taking into account these definitions we can state 
the following result of Walcher proved in [TU^ . 
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Theorem 15 |1U4] Let V{x,y) be a polynomial inverse integrating factor of a 
polynomial system x = P{x,y), y = Q{x,y) of degree d with P and Q coprime. 
Assume that the highest homogeneous components of P and Q of degree d are 
coprime. If there is an infinite simple critical point of the system, then the degree 
of V is exactly d + 1 . 

We observe that Theorem [15] is also proved in [12j with the additional assumption 
that the singularity at infinity p satisfies V{p) = where V is the projectivization 
of V. We remark that the approach introduced in |104j uses analytical techniques 
such as the Poincare-Dulac normal form and the proof given in [12] is completely 
algebraic and based on the extension of differential equations to the complex pro- 
jective plane and the results of Seidenberg about the reduction of singularities. 
The structure of polynomial inverse integrating factors is also studied by Walcher 
in [T05l[T06] . 

In relation with rational first integrals and in order to state the main result of 
[22] . we need to introduce some preliminary concepts, see also [38]. Let H = f /ghe 
a rational first integral of a polynomial system ([7]) . According to Poincare [92] we 
say that cGCU{cxo}isa remarkable value of H if f + eg is a reducible polynomial 
in C[x,?/] (here, c = oo means that f + eg denotes just g). In the work [22J it 
is proved that there are finitely many remarkable values for a given rational first 
integral H. 

Let now H he a polynomial first integral of degree n of a polynomial system ([7]). 
We say that the degree of H is minimal between all the degrees of the polynomial 
first integrals of ([7]) if any other polynomial first integral of ([7]) has degree > n. 

Assume H = f /g to he a rational first integral. Hence, we say that H has 
degree n if is the maximum of the degrees of / and g. Moreover, we say that 
the degree of H is minimal between all the degrees of the rational first integrals of 
system ([7]) if any other rational first integral of ([7]) has degree > n. 

Now suppose that c G C is a remarkable value of a rational first integral H = 
f / g and that YYi=i "^T factorization of the polynomial f + cg into irreducible 

factors in C[x, y]. If some of the is larger than 1, then we say that c is a critical 
remarkable value of H and that Ui = having > 1 is a critical remarkable 
invariant algebraic curve of ([7]) with exponent Oj. 

Finally, let / be a polynomial. We denote by / the homogeneous part of / of 
highest degree and this notation is also used for a Darboux functions like 

The main result of [22] is the following one. 

Theorem 16 [22] Suppose that a complex polynomial vector field X = P{x,y)dx + 
Q{x,y)dy of degree d with P and Q coprime has a Darboux first integral H given 
by where the polynomials fi and gi are irreducible and the polynomials gi and 
hi are coprime in C[x,?/]. Then the following statements hold. 
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(a) The inverse integrating factor V\ogH associated to the first integral log if is 
a rational function, and it can be written in the form ViogH = YliLi'^i' i^iU) 
with Ui G C[x, irreducible and ki G Z. Moreover, if X has no rational first 
integrals, then V\ogH is the unique rational inverse integrating factor of X . 

(b) Assume that H is a minimal polynomial first integral. Then there exists a 
polynomial inverse integrating factor. 

(c) Suppose that H = f/g is a minimal rational first integral of X and that X 
has no polynomial first integrals. It is not restrictive to assume that f and g 
are irreducible. Then, 

(c.l) the rational function 

V - 

where the product runs over all critical remarkable invariant algebraic 
curves Ui = having exponent ai is an inverse integrating factor; and 

(c.2) X has a polynomial inverse integrating factor if and only if H has at 
most two critical remarkable values. 

Additionally, if we assume for the first integral ([9]) that fi (respectively gj) are 
different for i = 1, . . . ,r (respectively j = !,...,£), and that it is complete (i.e. the 
unique algebraic invariant curves of system X are the fi = and the gj = 0), then 
the following two statements hold. 

(d) If X has no rational first integrals, then the inverse integrating factor V\ogH 
associated to the first integral log H is the polynomial 

V^ogH = flf.tlgr'- 

i=l j=l 

(e) If H is a multi-valued function and exp{hj / gj) are exponential factors of X 

for j = then ViogH = 111=1 /« 11^=1 d]^^^^ ^■^ ^ polynomial of degree 

d+ 1. 

In the particular case that /ij = for i = 1, . . . in the expression of ([9j), 
statement (d) of Theorem [16] can be thought as a generalization of following result 
due to Kooij and Christopher [77j and independently to Zolcidek ^109j . 

Theorem 17 [7711109] Consider a polynomial vector field X = P{x,y)dx+Q{x,y)d^ 
of degree d with P,Q E C[x,y] (resp. P,Q E M.[x,y]) having q invariant algebraic 
curves /« = such that the polynomials fi are irreducible in C[x,y] (resp. M.[x,y] 
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and satisfy that no more than two curves meet at any point of the plane (resp. 
M."^) and are not tangent at these points, no two curves have a common factor in 
their highest order terms and the sum of the degrees of the curves is d+ 1. Then, 
111=1 fi '^'^ inverse integrating factor of X and ni=i fl^ fof convenient Aj G C 
(resp. Xi E R) is a first integral of X . 

From statement (d) of Theorem (THl the following result easily follows. 

Corollary 18 [2^ Suppose that a real polynomial vector field X = P{x,y)dx + 
Q{x,y)dy of degree d with P and Q coprime has a Darhoux first integral H given 
by (Q where the polynomials fi and gj are irreducible, fi ^ gj, the polynomials 
gj and hj are coprime in R[x,y], exp{hj/gj^) are exponential factors of X , the Aj 
and Hj are either real numbers, or if some of them is complex then it appears its 

conjugate. If H is complete and H is multi-valued, then ViogH = 111=1 fi 11^=1 9]^~^^ 
is a polynomial of degree d + 1. // the system has foci or limit cycles, these are 
contained in the set {ViogH = 0}. 

Using Corollary [18] particularized to quadratic (d = 2) polynomial vector fields, 
in [22j it is obtained the next result. 

Corollary 19 [22j Under the assumptions of CorollarylTEl there are no real quadratic 
polynomial vector fields with a Darboux first integral ([9]) and a limit cycle. 

Some examples of polynomial systems satisfying the assumptions of Corollary 
[T8] are the following ones, see again [22] : 

• X = -y - xfi{x,y), y = X - yfi{x,y) where fi{x,y) = + - 1 = is 
an invariant circle which becomes an algebraic limit cycle. The origin is a 
focus and f2{x, y) = x"^ + y'^ is another invariant algebraic curve. The system 
possesses the inverse integrating factor V = fif2- 

• The system x = y — 4:xy, y = —x + + 2xy — y^ has the invariant algebraic 
curves fi{x,y) = 1 - 4x, f2{x,y) = V2y + {x + y - l)i and f?,{x,y) = 
\/2y — {x + y — l)i with = —1. The function V{x,y) = /1/2/3 is an 
inverse integrating factor. Notice that the system has a center at (0, 0) and 
a unstable focus at (1,0) 

The recent work [48j is also devoted to study the properties of remarkable val- 
ues. The polynomial R{x,y) := Hi'^r"^ defined as the product of all remarkable 
curves powered to their respective exponent minus one, is called the remarkable 
factor. From Theorem [121 if is a polynomial first integral, then the remarkable 
factor i? is a polynomial integrating factor of X. Moreover R divides the product 
n(if-l-Ci) where q are all the critical remarkable values of H. Thus the polynomial 



20 



Vji = Yli^ + Ci)/R is an inverse integrating factor of the system. The following 
theorem gives some relations between the degree of a system with a polynomial 
first integral, the degree of its inverse integrating factor Vr(x, y) and the number 
of critical remarkable values. 

Theorem 20 |18] Let X be a polynomial vector field of degree d and let H be a 
minimal polynomial first integral of X . Consider the remarkable factor R and the 
polynomial inverse integrating factor Vr. Let k be the number of critical remarkable 
values. Then, 

(a) k{k + d) < deg Vr = kdegH - degR < k{degH - 1) < degR{degR + d) 
and 

(b) deg Vr < deg H if and only if k = 1 . Moreover in this case deg Vr = d + 1. 



5 On the center problem 

One of the classical problems in the qualitative theory of planar analytic differential 
systems is to characterize the local phase portrait near an isolated singular point. 
By using the blow-up technique, see |13], this problem can be solved except when 
the singularity is monodromic, that is, it is either a focus or a center. The problem 
of distinguishing between a center or a focus is called the center problem. Another 
interesting problem is to know whether there exists or not a local analytic first 
integral defined in a neighborhood of a singular point. These two problems are 
equivalent when the singularity has associated nonvanishing complex conjugated 
eigenvalues. In this case, translating the singular point at the origin, after a linear 
change of variables and a rescaling of the time variable, the system can be written 
into the form: 

i = -y + f{x,y) , y = X + g{x,y) , (14) 

where f{x,y) and g{x,y) are analytic functions near the origin without constant 
nor linear terms. It is well known since Poincare and Liapunov that system (IT^ 
has a center at the origin if and only if there exists a local analytic first integral 
of the form H{x,y) = x"^ + y"^ + F{x,y) defined in a neighborhood of the ori- 
gin, where F starts with terms of order higher than 2. We recall here that the 
Poincare-Liapunov constants are the values V2k defined from the formal power 
series H{x,y) = Yl'^=2Hn{x,y), where H2{x,y) = (a;^ + ?/^)/2 and Hn{x,y) are 
homogeneous polynomials of degree n satisfying XH = YlT=2^'2k{x'^ + y^)''- The 
origin is a center of (HM if and only if all the Poincare-Liapunov constants vanish. 
When V2j = for j = 2, 3, . . . , /c — 1 and V2k 7^ 0, we say that the origin of system 
( IT4|) is a focus of order k. 
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The existence of invariant algebraic curves is strongly related with the origin 
of system flTll) being a center, as it is explained in [211 ESI llOU] . 



The proof of the following result is a particular case of a theorem that was given 
by Reeb in [95] (see also Mattel and Moussu |8l] and Moussu [87]). For a proof 
using elementary methods see [57] . 

Theorem 21 [95] System ( |T4l) /ias a center at the origin if and only if there is a 
nonzero analytic inverse integrating factor in a neighborhood of the origin. 

In fact, given a system (fT^ . the computational problems of looking for a first 
integral H{x, y) = x'^ + y'^ + - ■ ■ or for an inverse integrating factor V{x, y) = 1 + ■ ■ - , 
where the dots denote terms of higher order, are of the same difficulty. Thus, the 
inverse integrating factor offers an alternative to the solution of the center problem. 
In [21], it has been noticed that for many systems of type ((Hj) having a center at 
the origin there is an inverse integrating factor V with very simple properties which 
can be globally defined in all and which is usually a polynomial. By contrary, 
the first integral is, in general, a complicated expression that can not be written 
in terms of elementary functions. 

In particular, when in system f[T^ the functions / and g are both quadratic 
then there exists a polynomial inverse integrating factor of degree 3 or 5, see [TO] , 
whereas the first integrals are far more complicated, see [83]. When the functions 
/ and g in system (fT^ are both cubic homogeneous polynomials and the origin is 
a center, there exists a polynomial inverse integrating factor of degree at most 10, 
as it is also shown in [TU]. In [SS], the authors study cubic systems of the form (IT^ 
and give some sufficient conditions for the origin to be a center. This conditions 
come from the imposition to the system to have an inverse integrating factor. 

The work [30J is a survey on isochronous centers, that is, centers of the form 
f[T^ such that all the periodic orbits surrounding the origin have the same period. 
Many families of isochronous systems are listed and an explicit expression of an 
inverse integrating factor is given in each case. We include here a couple of results 
of the ones appearing in [30] which we have chosen for being the most known 
examples of isochronous centers in the literature. 

The quadratic systems with a isochronous center at the origin are characterized 
in the following result. 

Theorem 22 The origin is an isochronous center of a quadratic system (fT4|) if, 
and only if, the system can be brought by means of an affine change of coordinates 
and a rescaling of time, to one of the following four systems. For each case in the 
list we include the corresponding inverse integrating factor V{x,y). 

(Si) X = -y + x^-y^,y = x{l + 2y), with V{x, y) = (1 + 2yY . 
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{S2) X = -y + x^, y = x{l + y), with V{x,y) = (1 + yf ■ 

(53) X = -y- y = x{l- fy), with V{x,y) = (3 - 16y)(9 - 2% + 32x2). 

(5'4) X = -y + fx^ - y = x{l + fy), with V{x,y) = (3 + 8y){9 + 96y - 
256x2 + 128^2). 

Let us consider a cubic polynomial system of the form (H^ and let us assume that 
it contains no quadratic terms, that is, it is the sum of a linear system and a cubic 
homogeneous system. We say that such a system is cubic and with homogeneous 
nonlinearities. The following results characterizes which of these systems have an 
isochronous center at the origin. 

Theorem 23 The origin is an isochronous center of a cubic system with homoge- 
neous nonlinearities of the form (1141) if, and only if, the system can be brought by 
means of an affine change of coordinates and a rescaling of time, to one of the fol- 
lowing four systems. For each case in the list we include the corresponding inverse 
integrating factor V{x, y). 

{SD X = —y + x^ — Sxy"^, y = x + 3x'^y — y^, with V{x, y) = (x^ + y'^Y. 

(5*2) X = —y + x^ — xy"^, ij = X + x^y — y^, with V{x, ?/) = (1 + 2xyY . 

{SD X = -y + 3x^y, y = x - 2x^ + 9xy^, with V{x, y) = {1 - 3x^)*. 

{SD X = —y — 3x^y, y = x + 2x^ — Qxy"^ , with V{x, y) = {1 + 3x2)"^. 

In the works [101 [H] , Chavarriga writes system ([T^ in polar coordinates and 
studies the existence of inverse integrating factors polynomial in the radial vari- 
able. In [211 [251 126], the authors look for possible inverse integrating factors for 
polynomial vector fields of the form 

X = -yd^ + xdy + Xs , (15) 

where Xs is a polynomial homogeneous vector field of degree s > 2. In particular, 
they use the quasi-polar coordinates (i?, ip) where R = r*~^ and (r, ip) are the polar 
coordinates, that is, x = rcosy?, y = r sirup. Next, it is assumed the existence of 
an inverse integrating factor V{R, (p) of (|T5l) which is polynomial in the variable 
i?, that is, of the form V{R^(p) = J2^-QVi{{p)R^ with Voi^ip) = 1 and where Vi{(p) 
are homogeneous trigonometrical polynomials of degree i{s — 1). This assumption 
is clearly equivalent to impose an inverse integrating factor of the form V{R, ip) = 
nf=i(l + ^j(v^)-R)"S with ai G M. The authors try to solve the system of equations 
for the unknown functions Xi{ip) in the cases p = 1,2, 3. The case p = 1 is totally 
solved. If p = 2, only is solved the case ai =02 = (s + l)/(s — 1) ± 1/2 with 
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arbitrary s. Finally, when p = 3 the following two particular cases are investigated: 
«i = ct2 = «3 = 5/3 and either s = 2 or s = 3. 

In [27], some invariants are determined from which a formal first integral for 
system f|T5|) can be computed. Moreover, this technique is applied to the problem 
of determining the centers of polynomial vector fields f|T5|) . Recall that a complete 
classification of such centers is known when s = 2,3 but only partial results are 
known in the cases s = 4 and s = 5. 

Theorem [21] is used in |69j to find conditions to have a center. In this work, 
Gine proposes a formal power series V{x,y) = '^'^=oVn{x,y), where Vo{x,y) = 1 
and Vn{x,y) are homogeneous polynomials of degree n such that XV — VdivX = 
X]fcl2 '^2fc(a;^ + y^)'^, where X is given by (IT^ and the constants V2k are called the 
inverse integrating factor constants. Using the above mentioned quasi-polar co- 
ordinates {R, if) it is shown that, if the Poincare-Liapunov constants = for 
k = 1, . . . ,m and Vm+i 7^ 0, then Vm+i = — ((m + — 1) + 2)Vm+i- In this 
spirit, the paper |32| is concerned with the existence of a formal integrating factor 
of planar analytic system having a non degenerate focus or center at the origin and 
gives an algorithm to calculate the Poincare-Liapunov constants of any order. 

Given a real analytic planar vector field Xq with a center at pq G M^, in [62] 
the authors say that this center is limit of a linear type center if there exists a 
1-parameter family of analytic planar vector fields with e > 0, defined in a 
neighborhood of p^ G and having a non degenerate center at p^ for all e > 
sufficiently small. The main results of [62j are summarized as follows. 

Theorem 24 [02] Let Xq be a real analytic planar vector field with a center at 
Po eM.'^. Then, the following holds: 

(i) If po is a nilpotent center, then it is limit of a linear type center. 

(ii) // Pq is a Hamiltonian degenerate center, then it is limit of a linear type 
Hamiltonian center. 

(iii) Ifpo is a time-reversible degenerate center, then it is limit of a time-reversible 
linear type center. 

In the work [70j, Gine continues the study of the analytic centers which are 
limit of linear type centers. It is proved that if a degenerate center has an ana- 
lytic inverse integrating factor V{x,y) which does not vanish near the center, then 
this degenerate center is also the limit of a linear type center (changing the time 
variable). The idea is as follows. Assume V{x,y) is an analytic inverse integrat- 
ing factor of the analytic vector field Xq = P{x,y)dx + Q{x,y)dy such that (0,0) 
is a degenerate center and 1^(0,0) 7^ 0. Hence, the rescaled vector field Xq/V is 



24 



hamiltonian near the origin with analytic first integral H{x,y). Therefore, since 
the perturbed vector field A'e = Xq/V + e{—ydx + xdy) is Hamiltonian too, the 
origin becomes a linear type center of for all e 7^ 0. 

6 Limit cycles 

Let V{x,y) be an inverse integrating factor in the open set W C of a C^iU) 
planar vector field X. That is, the vector field X /V has zero divergence, where 
defined. If W is any simply connected component of W\\^~^(0), then the condition 
div(A'/^) = implies that X /V is Hamiltonian on W with single-valued 
hamiltonian function if : W — >^ M. Since Hamiltonian systems are area-preserving, 
hence have no limit cycles, and X and X /V are topologically equivalent, it follows 
immediately that, in the presence of an inverse integrating factor V , any limit 
cycle of X lies either in V~^{Q) or in a component ofU\ V^^{0) that is not simply 
connected. Using the machinery of de Rham cohomology, Giacomini, Llibre, and 
Viano eliminated the latter possibility in [53]. Hence, they prove the following 
theorem. 

Theorem 25 [63] Let j be a limit cycle of a real planar vector field X and let 
V be any inverse integrating factor of X defined in some neighborhood of'j. Then, 

7 c v-^{0). 

A different proof of Theorem [SS] can be found in [H \5U\ . We would like to recall 
here that M.V. Dolov in [315] studies the existence of a single valued regular inte- 
grating factor in a neighborhood of a limit cycle and presents some connections 
between an integrating factor and a limit cycle. Moreover, in the works [lOl BT] of 
Dolov and coauthors, published before the proof of Theorem [25l it is shown that 
vector fields with a Darboux inverse integrating factor of the form V = exp(i?) 
with rational R cannot have limit cycles. 

Theorem [25] has been applied in many papers to study limit cycles of a system 
as we will see in forthcoming sections. As an example where this theorem is ap- 
plied, we would like to recall the result of Llibre and Rodriguez in [81] where it is 
shown that every finite configuration of disjoint simple closed curves of the plane 
is topologically realizable as the set of limit cycles of a polynomial vector field. 
Moreover, the realization can be made by algebraic limit cycles, and an explicit 
polynomial vector field exhibiting any given finite configuration of limit cycles is 
provided. The proof of this realization makes use of the inverse integrating factor 
and, in particular, of Theorem [251 A generalization of the result of Llibre and 
Rodriguez is given in [HI] for systems in higher dimension, that is, it is shown that 
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any finite configuration of (smooth) cycles in R" can be realized (up to global dif- 
feomorphism) as hyperbolic and asymptotically stable limit cycles of a polynomial 
vector field. 

From Theorem [25] we have that the knowledge of an inverse integrating factor 
for a planar differential system ([3]) implies the knowledge of the number (and 
location) of the limit cycles of the system. Many authors have treated the problem 
of the existence of an inverse integrating factor. For a polynomial system (thus 
defined in the whole M^), the knowledge of a polynomial inverse integrating factor 
solves the question of the number and location of limit cycles of the polynomial 
system, see Section HI In the authors study the problem of existence of a 
polynomial inverse integrating factor in several cases of quadratic vector fields X. 
If such an integrating factor V{x, y) exists, then from Theorem [25] the curve V = 
is invariant for X and any limit cycle of X lies in this curve. Therefore, in [55] . 
the authors study planar quadratic polynomial vector fields that can have limit 
cycles and study the nonexistence of invariant algebraic curves, polynomial inverse 
integrating factors and algebraic limit cycles of arbitrary degree for these systems. 
Ye Yian-Qian |lU8j classified real quadratic systems that can have limit cycles in 
the following three families 

X = Sx — y + ix^ + Mxy + Ny^ , y = x{l + ax + by) , 

according to: family {!) if a = b = 0; family (//) if a 7^ and 6 = 0; family (///) 
if 6 7^ 0. In [55] it is proved that there are not algebraic limit cycles except for 
iN6 7^ and — 4iN > in family (/) (this result is improved in [15] where it 
is proved that there is no algebraic limit cycle for family (/)). Moreover, they also 
prove that the polynomial inverse integrating factors into families (/), (//)7v=0) 
(///)a=o and (//J)7v=o generically have at most degree 3. So, in the studied cases, 
the existence of polynomial inverse integrating factor implies the nonexistence of 
limit cycles or at most the existence of a circle as a unique limit cycle. 

Another interesting example of application of Theorem [2S] is given in the proof 
of several extensions to the Bendixson-Dulac Criterion to study of the number of 
limit cycles of planar differential systems, see [57] [58l [59] . An open set Z// C 
with smooth boundary is said to be ^-connected if its fundamental group, rriiU) 
is Z * . . * Z , OT in other words if U has i gaps. The classical Bendixson-Dulac 
Criterion is the following proposition, see [ST] for the statement and a short proof. 

Proposition 26 (Bendixson-Dulac Criterion) Let U he an open i-connected sub- 
set o/M^ with smooth boundary. Let X = P{x,y)dx + Q{x,y)dy be a vector field 
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of class C 



in lA. Let g : U ^ W be a C 



function such that 
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does not change sign in lA and vanishes only on a null measure Lebesgue set, such 
that {M = 0}r\{g = 0} does not contain periodic orbits of X . Then the maximum 
number of periodic orbits of X contained in U is I. Furthermore, each one of them 
is a hyperbolic limit cycle that does not cut {g = 0} and its stability is given by the 
sign of gM over it. 

7 The zero set of inverse integrating factors 

In Theorem [25] it is shown that hmit cycles 7 of a real planar vector field X 
belong to the zero set of any inverse integrating factor of X defined near 7, that 
is, 7 C V-^{0). 

In addition to containing any limit cycle of X lying in U, the zero set of V is also 
often connected to the separatrices of critical points of X in U. To understand why, 
recall that integral curves of X that map to themselves under the action of a Lie 
group are invariant solutions for the Lie group. Recall that when y = C,{x, y)dx + 
?7(x, y)dy is the infinitesimal generator of a nontrivial local Lie group of symmetries 
of X then the function V^(x, = det{A',3^} is an inverse integrating factor of X, 
as it has already been stated in Section O It is obvious that every solution of X 
which remains invariant under the action of the group with infinitesimal generator 
y must satisfy V{x, y) = 0. In other words, inverse integrating factors must vanish 
on invariant solutions. 

Based on these ideas, Bluman and Anco [B] argue heuristically that separatrices 
should also lie in V~^{0). Of course, any saddle loop in a Hamiltonian system is 
composed of separatrices not lying in the zero set of the trivial inverse integrating 
factor V = 1. Nevertheless, the idea has merit and we expect the zero set of V 
to play a role in the dynamics of X and it is very surprising that this fact was 
not completely accomplished until recent times. We repeat verbatim the following 
historical development on this issue given in [3J. 

"To our best knowledge, J. M. Page was the first author in making 
an observation of this kind. Concretely, the idea is developed in [88] of 
using Lie groups in the computation of singular solutions to the implicit 
first order differential equation F{x,y,y') = 0. The same idea is gath- 
ered in pgs. 113 and ss. of [88j, where several examples of calculation 
of envelopes are given, and later quoted without variations in [33], pgs. 
66 and ss.. In pg. Ill of his classical textbook [75], first published in 
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1926, E. L. Ince rescue Page's observation on envelopes but no other 
material is added. It would took several decades until some advance 
along this line of thinking might be registered. In this regard, the 
works of W. H. Steeb, C. E. Wulfman and G. D. Bluman and S. Kumei 
must be cited. In |102] . Steeb discussed the connection between limit 
cycles of two-dimensional systems and one-parameter groups of trans- 
formations. In |lU7j . Wulfman stated apparently general conditions on 
the infinitesimal generator of a Lie group admitted by a system of au- 
tonomous differential equations in order that an invariant solution is 
a limit cycle of the system. However, the argument he offers to sup- 
port these conditions rests on a heuristic more than rigorous basis. In 
turn, chapter 3 of contains a section (Section 3.6) devoted to dis- 
cuss the relationships existing between invariant solutions on one hand 
and exceptional paths on the other. Even though the developments in 
this section of the book seems to remain also on a semi- heuristic level, 
several examples and exercises are provided showing how the technique 
works in particular systems." 

In addition to this exhaustive historical description, we also would like to add 
the work of Gonzalez-Gascon [71j where it is pointed out that if there is an in- 
finitesimal generator of a Lie symmetry 3^ of a vector field X in W^, then on the 
limit cycles (periodic isolated orbits) of X it follows that X and y are parallel. 
This implies, in the particular case of planar fields that the associated inverse in- 
tegrating factor V = det{X, y} vanishes on the limit cycle. 

In [3] Berrone and Giacomini showed that, under mild additional hypotheses, 
the separatrices of hyperbolic saddle-points lying inU are contained in V^~^(0), and 
extended this result by showing that if F is a compact limit set all of whose critical 
points are hyperbolic saddle-points, then under mild conditions F C V~^{0) holds. 
Now, we summarize the results in [3]. 

It is easy to see that isolated vanishing points of an inverse integrating factor are 
singular points of the vector field. Moreover, for non-degenerate singularities (sin- 
gularities po € of A' with non vanishing Jacobian determinant det{DX{pQ)) ^ 0) 
one has the following result. 

Theorem 27 [3] Let po be a non-degenerate critical point of a vector field X 
and let V be an inverse integrating factor defined in a neighborhood of po and 
satisfying V{po) ^ 0. If det{DX{pQ)) > then po is a center. On the contrary, 
when det{DX{pQj) < 0, pq is a saddle-point. 

Next theorem is concerned with the stability of isolated zeroes of an inverse 
integrating factor. 
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Theorem 28 [3] Let po be an isolated zero of a non-negative inverse integrating 
factor V of a vector field X defined in a neighborhood U of pq. Then po is a 
stable (resp. unstable) singular point of X provided that divX\K < (resp. >). 
Furthermore, po is asymptotically stable (resp. unstable) provided that dwX\u < 
(resp. > 0). 

When a singularity po of A" is a non-isolated zero of an inverse integrating 
factor the following result holds. Here, given an orbit 70 of X, we denote by 0^(70) 
and a (70) its tu-limit set and a-limit set respectively. 

Theorem 29 [3] Let po be a non-isolated zero of an inverse integrating factor V 
of a vector field X . Then, one of the following two possibilities may occur: 

(a) There exists at least an orbit 70 of X (different of Pq) such that uj{'jo) = po 
or q;(7o) = Po and V\^g = 0. 

(b) There exists a infinite sequence {7n}neN of periodic orbits of X accumulating 
at Po such that V\-y^ = 0. 

A singularity Po of the vector field X is called strong if divA'(po) 7^ 0. Otherwise, 
when (HvX^pq) = 0, it is called weak. For a linear strong saddle points, it is easy to 
see that every inverse integrating factor must vanish on all four separatrix curves 
of the saddle. As it is established by the next theorem, the situation with nonlinear 
hyperbolic saddle points of systems is entirely analogous to the linear case. The 
proof is based on the normal form of X near the hyperbolic saddle po and the 
Stable Manifold Theorem. 

Theorem 30 [3] Let po be a hyperbolic saddle-point of a vector field X and V 
an inverse integrating factor defined in a neighborhoodU ofpQ. Then V vanishes on 
all four separatrix curves of the saddle provided that one of the following conditions 
holds: (i) Po is strong; (ii) po is weak and V{po) = 0. 

In the work [50], the previous theorem is slightly improved. If po is a hy- 
perbolic saddle point of a C''^^ vector field X whose k^^ saddle quantity is not 
zero and V is an inverse integrating factor defined in a neighborhood of po, then 
V{po) = (and, thus, V vanishes on all four separatrix curves of the saddle). For 
a full definition of saddle quantities see Subsection 18. 2l in relation with system fl23l) . 

As a corollary of Theorem [301 one can ensure the vanishing of an inverse inte- 
grating factor defined near certain saddle connections. Recall that a saddle con- 
nection is a union of saddle points and orbits connecting them. 
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Corollary 31 [5] Let V be an inverse integrating factor defined in a region con- 
taining a saddle connection T whose critical points are Pi for i = 1, . . . ,n. If V 
vanishes at a certain singular point pk, then V\t = 0. 

A graphic T = U*L]^0j(t) U {pi, . . . ,pk} is formed by k singular points pi, . . . ,Pk, 
Pk+i = Pi and k oriented regular orbits . . . , 0fc(t), connecting them such that 

4>i{t) is an unstable characteristic orbit of pi and a stable characteristic orbit of 
Pi+i- A graphic may or may not have associated a Poincare return map. In case it 
has one, it is called a polycycle. 

Now, let us suppose that F is a graphic, that is, F is a limit set which differs 
from a critical point or a periodic orbit. 

Theorem 32 |3] Let V be an inverse integrating factor defined in a region con- 
taining a compact graphic F. Then, the following holds: 

(a) V vanishes at a critical point at least ofT. 

(b) // all the critical points on F are non-degenerate, then V\r = 0. 

The main results of the paper [56j are generalizations and extensions of the 
previous results stated in [5]. A key ingredient in the proof of the results of ^B] is 
the concept of an integral invariant, introduced by Poincare in [93] for arbitrary di- 
mension, and its relation to inverse integrating factors. We denote by (l){t; {xo,yo)) 
the solution of ([3]) passing through the point {xo,yo) E U a.t t = 0; (j)(t; D) will 
denote the image of a domain D (ZU under the time-t map of the flow generated 
by the solutions of system 

Definition 33 Let fi : U G M."^ —>■ be a non-zero integrable function on U. The 
integral 

I y) dxdy (16) 

J<t>(t\D) 

is an integral invariant of system (I2D if for any measurable set D dU the integral 
is independent oft. 

The function /i is called the density of the integral invariant, based on the 
obvious hydrodynamic interpretation. Various versions of the following result can 
be found in textbooks, see for instance [2j. We state it in a form suited to our 
needs. In [56] is also provided a short proof. 

Lemma 34 [HS] Let U be an open subset o/M^, let V : U ^ be a function, 
and define a function fi : U \ ^^"'^(0) ^ M. by ^ = l/V . Then V is an inverse 
integrating factor of system in U if and only if the integral f|T6|) is an integral 
invariant for system onU \ V~^{0). 
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By using the relationship between inverse integrating factors and integral in- 
variants given in Lemma [Ml it is easy to see the next result. The definition of 
parabolic or elliptic sector can be found, for instance, in pHj . 

Theorem 35 |56] Let po be any critical point of system ([3]) at which there is an 
elliptic or parabolic sector. If V is any inverse integrating factor of ([3]) defined on 
a neighborhood of po, then V{po) = 0. 

In order to state the next result, we recall that a function / is called a Morse 
function if all its critical points are nondegenerate, i.e., the associated Hessian 
matrix has maximal rank at all the critical points. For Morse functions it is well 
known, see [71| for instance, that the set of critical points is discrete, that is, has 
no accumulation points. 

Theorem 36 Let T be any compact a- or u-limit set of system that 
contains a regular point, and let V be any inverse integrating factor of defined 
in some neighborhood of T. Depending on the smoothness of V , the following 
statements hold. 

(a) There exists a point p in T such that V{p) = 0. 

(b) // V is C"^ , then either T contains a point that is an accumulation point of 
isolated critical points ofV orTc V^^{0). 

(c) IfV is real analytic or Morse, then T C ^^"'^(0). 

Theorem [301 does not hold, in general, for non-hyperbolic singularities. But it 
does generalize for saddle or saddle-node singularities with exactly one non-zero 
eigenvalue, as the next result shows. 

Theorem 37 [56] Suppose po is an isolated singularity of a vector field X , and 
that V is an inverse integrating factor for X defined in a neighborhood of pq. If 
the linear part DX{pq) has exactly one zero eigenvalue, then V vanishes along any 
separatrix of X at Pq. 

We finish this section by stating a corollary of Theorems [301 and [371 

Corollary 38 [56] Let T be a polycycle (or graphic which need not be a limit 
set) of system ([3]) and let V be any inverse integrating factor of Q defined in 
some neighborhood of T . Assume that the critical points of ^ that belong to 
r are hyperbolic saddles pi,p2, . . . ,Pn or saddles and saddle-nodes qi,q2, ■ ■ ■ ,(lm 
with exactly one zero eigenvalue. If the separatrices of T are such that they always 
connect either p^ withpj and V{pk) = orp^ with qj or q^ with qj then T C V~^{0). 
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As the authors of remark, the hypothesis in Theorem [2S] that V be real 
analytic does not seem to be essential. Thus, in [SB] it is conjectured that only in the 
class for V, Theorem [36] remains valid. This conjecture was solved positively in 
^47j . In summary, in that paper the authors prove that there always exists a smooth 
inverse integrating factor in a neighborhood of a limit cycle and obtain a necessary 
and sufficient condition for the existence of an analytic one. This condition is 
expressed in terms of the Ecalle-Voronin modulus of the associated Poincare map. 
We recall that a germ of a map in the set of real analytic diffeomorphisms near the 
origin of M is analytically embeddable, i.e., it is the time-one map of an analytic 
vector field on the line, if and only if its Ecalle-Voronin modulus is trivial. The 
embedding properties of the Poincare map are crucial for the proof of the next 
theorem. 

Theorem 39 [17] Let 'j be a limit cycle of the analytic planar vector field X . Then 
there exists a neighborhood U of 'y and a function V G C^iU) which is an inverse 
integrating factor of X and vanishes exactly on 7. Moreover, V can be chosen 
analytic if and only if the Ecalle-Voronin modulus of the germ of the Poincare 
map of X along the limit cycle 7 is trivial. 

Corollary 40 [47J If '-y is a hyperbolic limit cycle of an analytic vector field X , 
then X admits an analytic inverse integrating factor in a neighborhood ofj. 

In addition, in ^47] it is also proved that a inverse integrating factor of a 
planar vector field must vanish identically on the polycycles which are limit sets of 
its fiow. We recall that a polycycle of a vector field is a compact invariant set 
which contains both regular and singular points. 

Theorem 41 [17] Let X be aC^ vector field defined in a domain U C M^. Suppose 
that r G U is a polycycle which is a limit set of X and X has a finite number of 
singular points in F. Then if X admits a inverse integrating factor V in U, 
then V C V~^{Q). 

The main idea of the proof of Theorem HT] is to pull back the vector field X and 
the inverse integrating factor V to the universal cover of U\{VX = 0} and exploit 
the fact that X /V lifts to a Hamiltonian vector field in the covering space. 

The existence of inverse integrating factors in a neighborhood of an elementary 
singularity is also established in [47j. The regularity of the inverse integrating 
factor depends on the kind of singularity and the proof makes crucial use of the 
theory of normal forms for planar vector fields. This considerably extends previous 
results in [2T] , where the authors prove for analytic vector fields the existence of a 
unique analytic inverse integrating factor in a neighborhood of a strong focus, or 
a non-resonant hyperbolic node, or a Siegel hyperbolic saddle. 
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The following result, which is stated in is a summary and a generalization 
of several results on the existence of a smooth and non-fiat inverse integrating 
factor Vo(x, y) in a neighborhood of an isolated singular point, see [211 El l67j . 

Theorem 42 Let the origin be an isolated singular point of (jS]) and let A, /i G C 
be the eigenvalues associated to the linear part of // A 7^ 0, then there exists 
a smooth and non-flat inverse integrating factor V{x, y) in a neighborhood of the 
origin. 

In |17] the existence of an analytic inverse integrating factor in a neighborhood 
of a non-degenerate monodromic singular point of an analytic system is charac- 
terized. If the origin is a non-degenerate center or a strong focus, there exists 
an analytic inverse integrating factor. If the origin is a weak focus, by Theorem 
|42] we have the existence of a smooth and non-flat inverse integrating factor, and 
there exists an analytic inverse integrating factor if and only if the Ecalle-Voronin 
modulus of the associated Poincare map is trivial, see also Theorem [391 In [47J the 
first known examples of real planar analytic vector fields not admitting an analytic 
inverse integrating factor in any neighborhood of either a limit cycle or an isolated 
singularity are given. 

In [51], we show the existence of an inverse integrating factor in a neighborhood 
of some degenerate singular points. 

Theorem 43 [5T] There exists an inverse integrating factor V^(x, y), of class at 
least , in a neighborhood of the following two types of singular points: a degen- 
erate focus without characteristic directions and a nilpotent focus. 

8 Bifurcations 

The inverse integrating factor has been shown to be very useful in many bifurcation 
problems. The books [73l [97] contain the main concepts and ideas of this theory 
in the framework of ordinary differential equations. 

Consider system ^ and take a parametric family of systems of the form 

X = V{x,y,e), y = Q{x,y,e), (17) 

where V{x, y, e) and Q{x, y, e) are analytic functions in (x, y) in the same open set 
as P{x, y) and Q{x, y) (or an open set we are interested in), are analytic for e near 
the origin and coincide with P{x,y) and Q{x,y) when e = 0, that is, V{x,y,0) = 
P{x,y) and Q{x,y,0) = Q{x,y). The parameter e is called bifurcation parameter 
and we assume it is defined in a neighborhood of the origin of M*^, with k e N; 
in many cases we consider that e is a real one-dimensional parameter {k = 1). 
For small values of the norm of £, we say that the family of systems (IT7I) is a 
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perturbation of system ([3]). When e takes values near the origin < << 1, the 
quahtative behavior of system (fT7|l can change with respect to the one of system ^ 
for e = 0. In this case, we say that a bifurcation has occurred. Bifurcation theory 
aims at characterizing under which conditions on system ([3]) and its perturbations, 
this bifurcations eventually happen and which are their properties. For example, 
consider a singular point p of system ^ and denote by A and /i the eigenvalues 
of the linearization of the system around p. If A, /i G M and A ■ yU < 0, then we 
say that p is a hyperbolic saddle and a classical result states that any perturbation 
of the system in a neighborhood of this point has the same qualitative behavior, 
that is, we have a saddle singular point that, when e tends to zero, tends to p. 
We give the adjective hyperbolic to those objects which maintain their qualitative 
nature under perturbations. In contrast, if p is a singular point of center type, i.e. 
it has a neighborhood filled with periodic orbits, then a perturbation of the system 
usually breaks these orbits and the point can be transformed, for instance, into a 
singular point of focus type, i.e. surrounded by orbits that spiral towards (or from) 
it. In this case, we say that p is a bifurcation point. When the considered family 
f[T7|) shows all the possible sample of qualitative behaviors that might occur when 
perturbing an object of system ([3]), we say that it is an unfolding. The minimum 
number of parameters needed to have an unfolding is called the codimension. 

Bifurcation theory is one of the most current tools used when trying to solve 
-^Qth YiiVoexi problem, part b. This problem was proposed in 1900 by D. Hilbert 
and asks for the maximum number and possible configurations of limit cycles that 
a polynomial system of the form ([7]) of degree d may have, only depending on the 
degree d. For a fixed system, Ecalle (1992) and Il'yashenko (1991) have demon- 
strated, in a different and independent way, that the number of limit cycles that 
the system may have is finite. However, the problem of determining whether there 
exists an upper bound on the number of limit cycles that a polynomial system of 
the form ([7]) can have, only depending on the degree d of the system, is still open. 

As R. Roussarie defines in [97], given a family of systems of the form (fT7I) . a 
limit periodic set is a compact and nonempty subset F of points so that there exists 
some succession (e„)„ which tends to when n +oo such that for every e„, the 
corresponding system (1171) has a limit cycle which tends to F, in the sense of 
the Haussdorf distance, when n +00. In this context, it is assumed that the 
parameters take values in a compact set. Following an analogous argument to the 
one used to prove Poincare-Bendixson Theorem, the structure of limit periodic 
sets can be determined. Given a limit periodic set of the family f|T7|) . we define 
its cyclicity as the maximum number of limit cycles which can be bifurcated from 
F in this family. In [HS], see also [97], R. Roussarie showed that the existence of 
a uniform upper bound in the number of limit cycles of an analytic family f[T7|) 
is equivalent to that each of its limit periodic sets F has finite cyclicity. This 
equivalence and the fact that all the limit periodic sets in a family (fT7I) can be 
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determined shows how bifurcation theory allows to tackle 16 Hilbert problem. In 
[46j, Dumortier, Roussarie and Rousseau established a list of 121 cases which are 
all the possible limit periodic sets that can appear within the family of quadratic 
systems and proposed a program, currently unfinished, to study all these graphics 
to demonstrate that there is a uniform upper bound for the number of limit cycles 
of polynomial systems of degree 2. 

The knowledge of inverse integrating factors for particular systems has simpli- 
fied its study and has allowed the understanding of several bifurcations. In [19] . 
for instance, the following family of cubic systems 

X = Xx — y + Xmix^ + (m2 — mi + mim2)x'^y + Xmim2xy'^ + m2y^, 
y = X + Xy — x"^ + Xmix^y + (mim2 — mi — l)xy'^ + Xmim2y^, 

where A, mi and m2 are arbitrary real parameters, is considered. The fact of 
knowing an inverse integrating factor 

V{x,y) := [x +y )(l + mia; + mim2y ) 

for this family of systems allows the determination of all the bifurcations within 
the family. 

Indeed, inverse integrating factors allow the understanding of the bifurcation of 
limit cycles from many limit periodic sets, as we explain in this section. The main 
result used in this context is Theorem [25] as it states that any inverse integrating 
factor defined in a neighborhood of a limit cycle needs to vanish on it. We recall 
that the zero set of a limit cycle is formed by orbits of the system. 

We split this section in three subsections depending on the considered limit 
periodic sets. 

8.1 Bifurcation from a period annulus 

In this subsection we consider planar differential systems of the form ^ with a 
singular point of center type. The set of periodic orbits surrounding this point is 
called its period annulus. A perturbation of the system usually breaks these peri- 
odic orbits but some of them might be maintained as limit cycles for the perturbed 
system. We say that this periodic orbits have bifurcated from the period annulus. 

There are several methods to determine how many limit cycles bifurcate from 
the periodic orbits of a period annulus. These methods are based upon differ- 
ent tools: the Poincare return map, see for instance p|; the Poincare-Pontrjagin- 
Melnikov integrals, see for instance [72]; averaging theory, see |98]; and the inverse 
integrating factor, see [Ml ESj, 1103] . This last method also gives the shape of the 
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bifurcated limit cycles up to any order of the perturbation parameter. 



We are going to explain the method described in J64j- Let us consider a Hamil- 
tonian planar differential system with a center at the origin: 

dH dH 

^ = iT^ y = ~ 

oy ox 

where H{x,y) is the Hamiltonian function and it is analytic in a neighborhood of 
the origin. We denote by V the period annulus of the center at the origin. Any 
analytic (nonzero) function of the Hamiltonian is an inverse integrating factor of 
the system. In particular, any (nonzero) constant function is an inverse integrating 
factor. 

Leu us consider an analytic perturbation of the previous Hamiltonian system: 
X = V{x,y,e), y = Q{x,y,e), (18) 

where 

V{x,y,e) ■■= ^ + fk{x,y), Q{x,y,e) := - ^ + J^^'' 9k{x,y), 

k=i k=i 

and where e is a small real parameter and fk{x,y), gk{x,y) are analytic functions 
in V U {(0, 0)}. Let us look for an analytic solution 



V{x,y,e) = ^e''Vk{x,y), 



k=0 

of the partial differential equation 

dx dy \dx dy J 

This partial differential equation gives a succession of linear differential equations 
for the functions Vk{x,y) which can be solved recursively. 

The equation of order in e implies that Vo{x,y) needs to be an inverse inte- 
grating factor for the unperturbed system. We have, thus, that Vq = Vo{H) is a 
function of the Hamiltonian H{x,y). 

The equation of order 1 in e gives a linear differential equation in Vi{x, y) whose 
non-homogeneous term contains the function Vo{x,y). Imposing that the function 
V needs to be periodic when evaluated on the unperturbed periodic orbits, it can 
be shown that Vo{h) needs to be 



Vo{h) = A / gi{x,y)dx - fi{x,y)dy, 

'{H=h} 
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where A is a nonzero real constant, and {H = h} denotes the periodic orbit of V 
contained in the h-level set of the Hamihonian H{x,y). Therefore, Vo(/i) is the 
first Poincare-Pontrjagin-Melnikov integral associated to system f|T8l) . Once we 
take this expression of Vq, we can solve the linear differential equation for Vi which 
is determined up to the sum of an arbitrary function of the Hamiltonian H. 

By induction on k it can be shown that Vk is determined up to the sum of an 
arbitrary function of the Hamiltonian H which we denote by Wk{h). When solving 
the linear differential equation of order k + 1 in e, and imposing that the function 
V needs to be periodic on the orbits of V, it can be shown that Wk{h) corresponds 
to the k + 1 Poincare-Pontrjagin-Melnikov integral associated to system (fT8|) . 

Indeed, in [M], the authors show that, fixed a small value of the zero sets 
of the functions '^^=0^'' ^k{x,y) give approximations up to order e"- of the limit 
cycles of system f|T8l) which bifurcate from V. When increasing the value of n, 
better approximations of these limit cycles are obtained and, thus, their shape is 
determined. 

In this method is generalized to non-Hamiltonian centers. The paper |103j 
purports a better understanding of this method as it studies this problem when 
the first i — 1 Poincare-Pontrjagin-Melnikov functions are identically zero. The 
main result in this paper is that, in this case, Vo{h) is the first non identically zero 
Poincare-Pontrjagin-Melnikov function. 

Most of these ideas are also used in [BB] to determine semistable limit cycles 
that bifurcate from V. Moreover, the method is applied to study the limit cycles 
which bifurcate from a Lienard system. 

We remark that this method is not only an alternative to the other methods 
as it shows how the inverse integrating factor is linked to bifurcation problems. 
This method is computationally as difficult as any other method but, moreover, it 
provides the shape of the bifurcated limit cycles. 

8.2 Bifurcation from monodromic cj-limit sets 

The work [50j is concerned with planar real analytic systems with an analytic 
inverse integrating factor defined in a neighborhood of a regular orbit (f){t). First 
of all it is shown that the inverse integrating factor defines an ordinary differential 
equation for the transition map along the orbit, see equation fl2Tl) . Taking two 
transversal sections Si and S2 based on (j)(t), it is studied the transition map of 
the flow of A' in a neighborhood of This transition map is studied by means 
of the Poincare map 11 : Ei — > E2. Given a point in Si, we consider the orbit of 
([3]) with it as initial point and we follow this orbit until it first intersects S2. 
Let {ip{s),ip{s)) G U, with s G X C M be a parameterization of the regular orbit 
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between the base points of Si and S2. Given a point {x,y) in a sufficiently 
small neighborhood of the orbit {ip{s),'ip{s)), we can always encounter values of 
the curvilinear coordinates {s, n) that realize the following change of variables: 
x(s,n) = ip{s) — mp'{s), y{s,n) = tp{s) + rnp'{s). We remark that the variable n 
measures the distance perpendicular to (f){t) from the point [x, y) and, therefore, 
n = corresponds to the considered regular orbit (f){t). We can assume, without 
loss of generality, that the transversal section Si corresponds to Si := {s = 0} 
and S2 to S2 := {s = L}, for a certain real number L > 0. We perform the change 
to curvilinear coordinates [x, y) 1— >■ (s, n) in a neighborhood of the regular orbit 
n = with s G X = [0,L]. Then, system (|3]) is written as the following ordinary 
differential equation: 

^ = F{s,n). (19) 
as 

We denote by \E'(s; uq) the flow associated to the equation f|T9l) with initial condition 
\E'(0; hq) = uq. In these coordinates, the Poincare map 11 : Si — > S2 between these 
two transversal sections is given by Il{nQ) = "^{L; Uq). 

We assume the existence of an analytic inverse integrating factor V{x, y) in 
a neighborhood of the considered regular orbit of the analytic system ([3]). 
In fact, when Si 7^ S2 and no return is involved, there always exists such an 
inverse integrating factor. The change to curvilinear coordinates gives us an inverse 
integrating factor for equation (HM . denoted by V{s,n) and which satisfies 

dV dV dF ~ 

^ + ^F{s,n) = ^V{s,n). (20) 
OS on on 

Now, we can state one of the main results of [50]. 

Theorem 44 |50j We consider a regular orbit (f){t) of the analytic system ([3]) which 
has an inverse integrating factor V{x,y) of class defined in a neighborhood 
of it and we consider the Poincare map associated to the regular orbit between 
two transversal sections 11 : Si — > S2. We perform the change to curvilinear 
coordinates and we consider the ordinary differential equation f ll9l) with the inverse 
integrating factor V {s , n) which is obtained from V{x,y). In these coordinates, the 
transversal sections can be taken such that Si := {s = 0} and S2 := {s = L}, for 
a certain real value L > 0. We parameterize T^i by the real value of the coordinate 
n. The following identity holds. 

y(L,n(n)) = •i/(0,n)n'(n). (21) 



Theorem HH is the key point to prove Theorems HH] and | 
Further, in [50j the authors consider regular orbits whose Poincare map is a 
return map and take profit from the result stated in Theorem HH in order to study 
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the Poincare map associated to a limit cycle or to a homoclinic loop, in terms of 
the inverse integrating factor. To do that, the following definition of vanishing 
multiplicity of an analytic inverse integrating factor V{x, y) of the analytic system 
([3]) over a regular orbit (f){t) is needed. 

Definition 45 [50j Let V{x,y) be an analytic inverse integrating factor of the 
analytic system ([3]) and (pit) a regular orbit of it parameterized by {ip{s),'ip{s)) e U, 
with s G X C M. Consider the local change of coordinates x{s, n) = (p{s) — mp'{s), 
y{s,n) = ip{s) + rnp'{s) defined in a neighborhood of the considered regular orbit 
n = and take the following Taylor development around n = 0: 

V{x{s,n),y{s,n)) = n'^vis) + 0(n™+^), (22) 

where m is an integer with m > and the function v{s) is not identically null, we 
say that V has multiplicity m on (f)(t). 

In fact, in [50] it is proved that v{s) ^ for any s G X, and thus, the vanishing 
multiplicity of V on (pit) is well-defined over all its points. 

Let us consider as regular orbit a limit cycle 7 and we use the parameterization 
of 7 in curvilinear coordinates (s,n) with s G [0,L). Thus, the Poincare map 
associated to 7 is Il{nQ) = \1/(L;?t,o). It is well known that 11 is analytic in a 
neighborhood of = 0. We recall that the periodic orbit 7 is a limit cycle if, and 
only if, the Poincare return map 11 is not the identity. If 11 is the identity, we have 
that 7 belongs to a period annulus. We recall the definition of multiplicity of a 
limit cycle: 7 is said to be a limit cycle of multiplicity 1 if n'(0) 7^ 1 and 7 is said to 
be a limit cycle of multiplicity m with m > 2 if n(no) = uq + jSmU^ + 0{n^^^) 
with f3m 7^ 0. Then, one has the following result for limit cycles. 

Theorem 46 [50] Let 'j be a periodic orbit of the analytic system ([3]) and let V be 

an analytic inverse integrating factor defined in a neighborhood of'j. 

(a) If '-f is a limit cycle of multiplicity m, then V has vanishing multiplicity m 
on 7. 

(b) IfV has vanishing multiplicity m on 7, then is a limit cycle of multiplicity 
m or it belongs to a continuum of periodic orbits. 

Since the Poincare map of a periodic orbit is an analytic function and the 
multiplicity of a limit cycle is a natural number, the following corollary is obtained. 

Corollary 47 [50j Let j be a periodic orbit of the analytic system ([3]) and let V be 
an inverse integrating factor of class defined in a neighborhood of'-/. We take the 
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change to curvilinear coordinates x{s,n) = ip{s) — nip\s), y{s,n) = ip{s) + rnp'{s) 
defined in a neighborhood of If we have that the leading term in the following 
development around n = 0: 

V{x{s,n),y{s,n)) = v{s) + o(n''), 

where v{s) ^0 is such that either p = or p > 1 and p is not a natural number, 
then 7 belongs to a continuum of periodic orbits. 

A regular orbit = {x{t),y{t)) of is called a homoclinic orbit if 0(t) po 
as t — s> ±oo for some singular point po. A homoclinic loop is the union F = 
(pit) U {po}- We assume that po is a hyperbolic saddle, that is, a critical point of 
system ([3]) such that the eigenvalues of the Jacobian matrix DX{pq) are both real, 
different from zero and of contrary sign. We remark that this type of graphics 
always has associated (maybe only its inner or outer neighborhood) a Poincare 
return map 11 : S ^ S with E any local transversal section through a regular 
point of r. We will assume that F is a compact invariant set. A goal in [50] is to 
study the cyclicity of the described homoclinic loop F in terms of the vanishing 
multiplicity of an inverse integrating factor. Roughly speaking, the cyclicity of F 
is the maximum number of limit cycles which bifurcate from it under an analytic 
perturbation of the analytic system ([3]). Before state the result for homoclinic 
loops, we recall briefly that the first saddle quantity is ai = divX{po) and it 
classifies the point po between being strong (when ai ^ 0) or weak (when ai = 0). 
If Po is a weak saddle point, the saddle quantities are the obstructions for it to 
be analytically orbitally linearizable. In order to define the next saddle quantities 
associated topo, we translate the saddle-point po to the origin of coordinates and we 
make a linear change of variables so that its unstable (resp. stable) separatrix has 
the horizontal (resp. vertical) direction at the origin. Let po be a weak hyperbolic 
saddle point situated at the origin of coordinates and whose associated eigenvalues 
are taken to be ±1 by a rescaling of time, if necessary. Then, it is well known the 
existence of an analytic near-identity change of coordinates that brings the system 
into: 

k-l 

X = X + ^ flj + Ok x^^^y^ + • • • , 

(23) 

y = -y - 5]aixV+^ - fofea;'^?/'^-*-^ + ■■■ , 

i=l 

with ttk — bk 7^ and where the dots denote terms of higher order. The first 
non- vanishing saddle quantity is defined by ak+i := au — bk- 

Theorem 48 [50] Let T be a compact homoclinic loop through the hyperbolic saddle 
Po of the analytic system ([3]) whose Poincare return map is not the identity. Let 
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V be an analytic inverse integrating factor defined in a neighborhood of T with 
vanishing multiplicity m over T. Then, m > 1 and the first possible non-vanishing 
saddle quantity is am- Moreover, 

(i) the cyclicity ofT is 2m — 1, if am 7^ 0, 

(ii) the cyclicity of T is 2m, otherwise. 

In addition, in [30] it is described one obstruction to the existence of an analytic 
inverse integrating factor defined in a neighborhood of certain homochnic loops. 
First of all, we recall some concepts. By an affine change of coordinates, in a 
neighborhood of a hyperbolic saddle, any analytic system can be written as i; = 
Xx + f{x,y), y = fiy + g{x,y), where / and g are analytic in a neighborhood of the 
origin with lowest terms at least of second order and /i < < A. This hyperbolic 
saddle is analytically orbitally linearizable if there exists an analytic near-identity 
change of coordinates transforming the system to a; = Xxh{x,y), y = fiyh{x,y) 
with h{0, 0) = 1. On the other hand, when /i/A = —q/p G Q~ with p and q natural 
and coprime numbers, the saddle is called p : q resonant. 

Proposition 49 [50] Suppose that the analytic system has a homoclinic loop 
r through the hyperbolic saddle point po which is not orbitally linearizable, p : q 
resonant and strong {p q)- Then, there is no analytic inverse integrating factor 
V{x,y) defined in a neighborhood ofT. 

8.3 Generalized Hopf Bifurcation 

Let us consider a planar real system ([3]), x = P{x,y), y = Q{x,y) and suppose 
that it is analytic near an isolated monodromic singular point po which we assume 
to be at the origin. We associate to system ([3]) the vector field Xq = P{x,y)dx + 
Q{x,y)dy. We consider an analytic perturbation of system ([3]) of the form: 

X = P{x,y) + P{x,y,e), y = Q{x,y) + Q{x,y,e), (24) 

where e G is the perturbation parameter, < ||£:|| << 1 and the functions 
P{x,y,e) and Q{x,y,e) are analytic for {x,y) G U, analytic in a neighborhood of 
e = and P{x, y, 0) = Q{x, y, 0) = 0. We associate to this perturbed system 
the vector field a; = {P{x,y) + P{x,y,€))d^ + iQ{x,y) + Q{x,y,e))dy. 

We say that a limit cycle 7^ of system (12^ bifurcates from the origin if it tends 
to the origin (in the Hausdorff distance) as e 0. We are interested in giving 
a sharp upper bound for the number of limit cycles which can bifurcate from the 
origin po of system under any analytic perturbation with a finite number p of 
parameters. The word sharp means that there exists a system of the form 
with exactly that number of limit cycles bifurcating from the origin, that is, the 
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upper bound is realizable. This sharp upper bound is called the cyclicity of the 
origin po of system and will be denoted by Cycl(A'£,po) all along this section. 

In |5T] we consider systems of the form ([3]) where the origin po is a focus 
singular point of the following three types: non-degenerate, degenerate without 
characteristic directions and nilpotent. The results of [51] do not establish that 
the cyclicity of this type of singular points is finite but give an effective procedure to 
study it. In the three mentioned types of focus points, we will consider a change to 
(generalized) polar coordinates which embed the neighborhood U of the origin into 
a cylinder C = {(r, ^) G M x 5^ : |r| < 6} for a certain sufficiently small value of 
6 > 0. This change to polar coordinates is a diffeomorphism in W — {(0,0)} and 
transforms the origin of coordinates to the circle of equation r = 0. In these new 
coordinates, system can be seen as a differential equation over the cylinder C 
of the form: 

To = ^C-")- P^' 

where J-'{r, 9) is an analytic function in C . We consider an inverse integrating 
factor V{r, 6) of equation (!25l) . that is, a function y : C — > M of class C^{C), which 
is non locally null and which satisfies the following partial differential equation: 

We remark that since V{r, 9) is a function defined over the cylinder C it needs to 
be T-periodic in 6', where T is the minimal positive period of the variable 9, that 
is, we consider the circle 5^ = M/[0,T]. The function V{r,9) is smooth (C°°) and 
non-fiat in r in a neighborhood of r = 0. 

Let us consider the Taylor expansion of the function V{r,9) around r = 0: 
V{r,9) = Vm{9)r"^ + 0{r"^~^^), where Vm{9) ^ for 9 E and m is an integer 
number with m > 0. We say that m is the vanishing multiplicity of V{r, 9) on 
r = 0. The uniqueness of V{r, 9) implies that the number m corresponding to the 
vanishing multiplicity of V{r, ^) on r = is well-defined. 

We consider a system of the form: 

X = Pdix,y) + Pd+iix,y), y = qdix,y) + Qd+i{x,y), (26) 

where (i > 1 is an odd number, pd{x,y) and qd{x,y) are homogeneous polynomials 
of degree d and Pd+i{x,y),Qd+iix,y) e C(||(a:, i/)!!*^"^^). We assume that p'^{x,y) + 
q'^(x,y) ^ 0. A characteristic direction for the origin of system fl26|) is a linear 
factor in M[a;, of the homogeneous polynomial xqd{x,y) — ypd{x,y). If there 
are no characteristic directions, then the origin is a monodromic singular point of 
system fl26|) . that is, it is either a center or a focus. 
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In relation with system (^^, an analytic perturbation field y, e), Q{x, y, e)) 
is said to have subdegree s if {P{x,y,e),Q{x,y,e)) = 0{\\{x,y)\\'^). In this case, 
we denote by A'i^' = {P{x, y) + P{x, y, e))dx + {Q{x, y) + Q{x, y, s))dy the vector 
field associated to such a perturbation. 

Theorem 50 [51] We assume that the origin po of system fl26l) is monodromic and 
without characteristic directions. Take polar coordinates x = r cos 6, y = r sin 6 
and let V{r, 6) he an inverse integrating factor of the corresponding equation fl2^ 
which has a Laurent expansion in a neighborhood of r = of the form V{r,9) = 
VmiO) r™ + 0(r™+i), with v^{e) ^ and m e Z. 

(i) Ifm<Oormis even, then the origin of system f l26p is a center. 

(ii) // the origin of system f l26p is a focus, then m > 1, m is an odd number and 
the cyclicity Cyd{X^,pQ) of the origin of system ( B^ satisfies Cyc\{X^,Po) > 
{m + d)/2 — l. In this casern is the vanishing multiplicity ofV{r,6) onr = 0. 

(11.1) //, moreover, the focus is non-degenerate {d = 1), then the aforemen- 
tioned lower bound is sharp, that is, CycI{X^,po) = (m — l)/2. 

(11.2) // only perturbations whose subdegree is greater than or equal to d are 
considered, then the maximum number of limit cycles which bifurcate 
from the origin is (m — l)/2, that is, Cycl(A:'i°'',po) = {rn — l)/2. 

Remark 51 From the proof of Tfeeorem lSOj, it follows that if there exists an inverse 
integrating factor Vo{x , y) of system (!26|) such that Vo(r cos6', rsin^)/r'^ has a Lau- 
rent expansion in a neighborhood ofr = 0, then the exponents of the leading terms 
of VqIt cos 6, rsm6)/r'^ and V{r,6) coincide. Thus, the vanishing multiplicity m 
can be computed without passing the system to polar coordinates. 

We assume that the origin of system (!26l) is a focus without characteristic 
directions and that the vanishing multiplicity of an inverse integrating factor on it 
is m. If system is written as i; = P{x,y) and y = Q{x,y), then the system: 

X = P{x,y) + xK{x,y,e), y = Q{x,y) + yK{x,y,e), (27) 

fc-i 

where L := {m -\- d)/2 — 1 and K{x,y,e) = e:*^"* (x^ + y^)*"' has at least 

(m + c/) /2 — 1 limit cycles bifurcating from the origin for convenient values of the 
real parameters Oj. We recall that both m and d are odd and d > 1, m > 1. 

We say that the origin of system ([3]) is a nilpotent singular point if it is a 
degenerate singularity that can be written as: 

X = y + P2{x,y) , y = Q2ix,y) , (28) 
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with P2{x,y) and Q2{x,y) analytic functions near the origin without constant 
nor hnear terms. The following theorem is due to Andreev [T] and it solves the 
monodromy problem for the origin of system fl28l) . 

Theorem 52 |T] Let y = F{x) be the solution ofy + P2{x, y) = passing through 
(0, 0). Define the functions f{x) = Q2{x, F{x)) = ax" + ■ ■ ■ with a ^ and a > 2 
and = {dP2/dx + dQ2/dy){x, F{x)) . We have that either = + • ■ • 
with b ^ and (3 > 1 or (f){x) = 0. Then, the origin of (!28|) is monodromic if, and 
only if, a < 0, a = 2n — 1 is an odd integer and one of the following conditions 
holds: (i) (3 > n — 1; (ii) (3 = n — 1 and 6^ + Aan < 0; (Hi) (f){x) = 0. 



Definition 53 We consider a system of the form fl28|) with the origin as a mon- 
odromic singular point. We define its Andreev number n > 2 as the corresponding 
integer value given in Theorem [52l 

We consider system (1281) and we assume that the origin is a nilpotent mon- 
odromic singular point with Andreev number n. Then, the change of variables 

{x,y) ^ {x,y-Fix)), (29) 

where F{x) is defined in Theorem |52l and the scaling 

(x,y) ^ i^x,-^y), (30) 

with ^ = (— l/a)^/*^^~^"\ brings system fl28|) into the following analytic form for 
monodromic nilpotent singularities 

X = y{-l+Xi{x,y)), y = f{x) +y(f){x) +y'^Yo{x,y), (31) 

where Xi(0,0) = 0, f{x) = x'^"'~'^ + ■■■ with n > 2 and either 0(x) = or 
4>{x) = bx^ + ■ ■ ■ with (3 > n — 1. We remark that we have relabelled the functions 
/(x), (j){x) and the constant b with respect to the ones corresponding to system 



We are going to transform system fl?Il) to an equation over a cylinder of the 
form (l25|l . The transformation depends on the Andreev number n and it is given 
through the generalized trigonometric functions defined by Lyapunov as the unique 
solution x{9) = Cs9 and y{9) = Sn0 of the following Cauchy problem 

^ = -y, I = ^^'^-S x(0) = 1, y(0) = 0. (32) 

We introduce in ]R^\{(0,0)} the change to generalized polar coordinates, {x,y) ^ 
(r, 9), defined by 

X = rCs9, y = r'^Sn^. (33) 

We consider the following definition, which will be used in the following Theo- 
rem! 
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Definition 54 An analytic perturbation vector field {P{x,y,e),Q{x,y,e)) is said 
to be (1, n)-quasihomogeneous of weighted subdegrees {wx,Wy) if P{\x, X^y,e) = 
OiX""^) and_Q{\x,X'y,e) = C(A^^). In this case, we denote by A']"'"-"'''! = 
{P{x,y) + P{x,y,e))dx + {Q{x,y) + Q{x,y,e))dy the vector field associated to 
such a perturbation. 

The following theorem is one of the main results of [ST] . The symbol [xj denotes 
the integer part of x. 

Theorem 55 |5T] We assume that the origin of system fl28l) is monodromic with 
Andreev number n. Take generalized polar coordinates (3^) and let V{r, 6) be an 
inverse integrating factor of the corresponding equation ( l25l) which has a Laurent 
expansion in a neighborhood of r = of the form V{r,6) = Vm{0)r"^ + 0{r"^~^^), 
with Vm{0) ^ and m G Z. 

(i) If m < or m + n is odd, then the origin of system is a center. 

(ii) // the origin of system fl25]) is a focus, then m > 1, m + n is even and its 
cyclicity Cycl{Xi;,po) satisfies Cjc[{Xs,po) > (m + n)/2 — 1. In this case, 
m is the vanishing multiplicity ofV{r,6) onr = 0. 

(iii) // the origin of system fl3Tl) is a focus and if only analytic perturbations of 
{1, n) -quasihomogeneous weighted subdegrees {wx,Wy) with > n and Wy > 
2n — 1 are taken into account, then the maximum number of limit cycles 



which bifurcate from the origin is [(m — l)/2j, that is, Cjd^Xe' " ,Po) = 



Remark 56 The proof of this theorem shows that if there exists an inverse in- 
tegrating factor VQ{x,y) of system (131 p such that VQ*(rCs^, r"'Sn6')/r^"'~^ has a 
Laurent expansion in a neighborhood of r = 0, then the exponents of the lead- 
ing terms ofV^{rCs6, r"" Sn 6) / r'^'"'^^ and V{r,6) coincide. Therefore, the value 
of m can be determined without performing the transformation of the system to 
generalized polar coordinates. 

We assume that the origin of system (13T|) is a focus with Andreev number n 
and that the vanishing multiplicity of an inverse integrating factor on it is m. If 
system (13T|) is written as x = P{x,y) and y = Q{x,y), then the system: 




L(m-1)/2J. 



X = P{x,y) + xK{x,y,e), 



y 



Q{x,y) +nyK{x,y,e), 



(34) 



where 



L-l 



K{x,y,e) = ^ 



€ 



L-i 



,2i 



i=0 
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and L = (m + ?2)/2 — 1, has at least (m + n)/2 — 1 limit cycles bifurcating from 
the origin for convenient values of the real parameters Oj. We recall that m and n 
have the same parity. 

The following corollary establishes a necessary condition for system fl28l) to 
have an analytic inverse integrating factor Vo(a;,?/) defined in a neighborhood of 
the origin. 

Corollary 57 We assume that the origin of system (!28|) is a nilpotent focus with 
Andreev number n, and that there exists an inverse integrating factor Vo(x, y) of 
(!28|) which is analytic in a neighborhood of the origin. Then, n is odd. 

8.4 Singular perturbations 

As some recent research papers show, see |15], limit periodic sets containing an 
infinite number of critical points may have a cyclicity higher than expected. Due 
to the narrow relationship between limit cycles and the inverse integrating factor, 
the context of singular perturbations is a brand new and very interesting place to 
apply properties of the inverse integrating factor in order to detect limit cycles. As 
we have seen, for other limit periodic sets, the inverse integrating factor does not 
only give an alternative way to study the cyclicity but contains more information: 
location of limit cycles, direct computation of the cyclicity of the object under 
study, an explicit partial differential equation (jl]) which gathers all the informa- 
tion, . . . 

As far as we know, the only work where the inverse integrating factor is related 
with a singular perturbation problem is [72j , where one-parameter families of vector 
fields A'e in of the form = f{x, y, e) dx + £g{x, y, e) dy, where e >0 and /, g 
are analytic functions, are taken into account. The aim of the singular perturbation 
problems is to study the phase portrait, for e sufficiently small, near the set of 
singular points of Xq, that is, S = {{x,y) G : f{x,y,0) = 0}. In particular, 
the question is to decide if X^ has a limit cycle which tends to a singular orbit of 
Xq when e \ 0. A singular orbit (also denoted as slow- fast cycle, see |l6j) is 
a limit periodic set of the system Xq. For the vector field Xq, we say that a point 
n G S is normally hyperbolic if {df /dx){n,0) ^ 0. The system of differential 
equations associated to X^ is 

X = f{x,y,e), y = eg{x,y,e), 

where the dot denotes derivation with respect to the time t. We call this system 
the fast system. By the time rescaling r = et, we get the slow system: 

ex' = f{x,y,e), y' = g{x,y,e), 
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where ' denotes derivation with respect to r. The reduced problem is defined by 
the slow system taking e = 0, which gives one differential equation constrained to 
the slow manifold or critical curve S, that is, the reduced problem is 

f{x,y,0) = 0, y' = g{x,y,e). 

The only singular orbits taken into account in [79j are the ones consisting of three 
pieces of smooth curves; an orbit of the reduced problem starting at a normally 
hyperbolic point ni G S, an orbit of the reduced problem ending at a normally 
hyperbolic point n2 G S and an orbit of the fast problem connecting the two 
previous ones. The main results are the following. 

Theorem 58 [79] Consider Eq > and an inverse integrating factor of 

Xg, that is Xi;{Vs{x,y)) = divA^ ^(x, y), defined in an open set U C for any 
< e < ^o- Let T G U be a singular orbit and be a limit cycle of X^, in U for 
e G (0,£:o); with T , according to the Hausdorff distance. Then Vq{V) = 0. 

Corollary 59 [79] Consider Ve{x,y) an inverse integrating factor of as in The- 
orem [581 If the level zero of the function Vo{x,y) does not contain a closed curve, 
then there exists > such that does not present a limit cycle for < e < eo 
in U. 

As an application of these results, the following examples are given in [79] . The 
following vector fields present no limit cycles because the corresponding inverse 
integrating factors have no closed curves in their level zero sets. 

• The vector field Xf, = {y"^ — x'^)dx + ex"^ dy has the inverse integrating factor 

V^(x, y) = y^ — yx"^ — x^e. 

• The vector field A'^ = {y — x'^)dx + e x dy has the inverse integrating factor 

V,{x,y) = -y + x' + (1/2)6. 

• The vector field X^ = {~y + x'^)dx + exdy has the inverse integrating factor 

V,{x,y)=y-x^ + (1/2)6. 

9 Some generalizations 

9.1 The inverse Jacobi multiplier 

Inverse Jacobi multipliers are a natural generalization of inverse integrating factors 
to n-dimensional dynamical systems with n >3. In [1], it is developed the theory 
of inverse Jacobi multiplier from its beginning in the formal methods of integration 
of ordinary differential equations to modern applications. 
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In this section we will assume that X = ^"^^^ Xi{xi, . . . , Xn)dxi is a vector 
field defined in the open set C R"-. A function : W — M is said to be 
an inverse Jacobi multiplier for the vector field X mU when V solves in U the 
linear first order partial differential equation XV = VdivX. The first appearance 
of these multipliers occurs in the works of C.G.J. Jacobi, about the middle of the 
past century. Many properties of inverse integrating factors for the planar case 
{n = 2) are inherited by inverse Jacobi multiplier. We list some of them: 

• If the change of coordinates y = (f){x) is introduced, then W{y) = (V o 
(f)^^){y) det{D(j){(j)^^{y))} is an inverse multiplier of the transformed vector 
field (l)^:X. 

• Let Vi and V2 be two linearly independent inverse Jacobi multipliers of X 
defined in U. If Vi(x) 7^ for all x G W, then the ratio V2/V1 is a first integral 
of X in U. 

• One can use local Lie groups of transformations to find inverse Jacobi mul- 
tipliers as follows. Assume X admits in W a (ra — l)-parameter local Lie 
group of transformations with infinitesimal generators {3^i, . . . , 3^„_i}. Then, 
an inverse Jacobi multiplier V for in is furnished by the determinant 



• Let {3^1, . . . , 3^n-i} be the generators of n — 1 local Lie groups of symmetries 
admitted by X in U. Then, the inverse multiplier V = det{X, yi, . . . , 3^n-i} 
vanishes on every invariant solution of X contained in U. Recall here that 
an invariant solution of X corresponding to the group G is defined to be an 
integral curve of X which is invariant under the action of G. 

• Let pq & U he an isolated zero of an inverse Jacobi multiplier V such that 
y > in a neighborhood Af of po. Then, pq is a stable (resp. unstable) 
singular point of X provided that divX < (resp. > 0) in Af. Further- 
more, the stability (resp. unstableness) of po is asymptotic stability (resp. 
unstableness) provided that divA" < (resp. > 0) in Af. 

In the following, we summarize some of the results obtained in [1]. By a limit 
cycle 7 of A* we mean a T-periodic orbit which is a or u;-limit set of another orbit 
of X. Let V be an inverse Jacobi multiplier defined in a region containing 7. If 
7 = {7(t) eU -.0 <t <T},we define 



As it is well known, A (7) is the sum of the characteristic exponents of the limit 
cycle 7. We recall that if A (7) > then 7 is not orbitally stable. We will say that 



v = det{x,yi,...,yr,.i}. 
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7 is a strong limit cycle when A(7) 7^ 0. If, on the contrary, A(7) = 0, then we 
say that 7 is a weak limit cycle. 



Theorem 60 [Ij Let V be an inverse Jacobi multiplier of X defined in a region 
containing a limit cycle 7 of X . Then, 7 is contained in l^~^(0) in the following 
cases: (i) if 'j is a strong limit cycle, or (ii) if 7 is asymptotically orbitally stable 
(unstable). 

Theorem 61 |1] Let V be a Jacobi inverse multiplier defined in a neighborhood of 
a limit cycle 7 of the vector field X. If '-j is a strong limit cycle, then 

(i) V vanishes on 1^(7), the stable manifold of'-f, provided that A(7) > 0; 

(ii) V vanishes on W'^{'y), the unstable manifold of'y, provided that A(7) < 0. 



The following example appears in fl]. Consider the cubic polynomial vector 
field in 



where f{x,y) = \ — — y^ and A > is a real parameter. The circle 7 = 
{/(x, ?/) = 0} n {2; = 0} is a limit cycle of system (!35l) of period T = 2n/\. In fact, 
7(t) = (cos At, sin At, 0). It is easy to compute that 



Jo ^ 
An inverse Jacobi multiplier for this system is 

Vi{x,y,z) = f{x,y){x'^ + y'^)z . 

In addition, when A = —1/2, then V2{x, y, z) = (x^ +1/^)^ is another inverse Jacobi 
last multiplier for system (1351) . 

As usual, a hyperbolic singular point po of a vector field X is named a saddle 
point when the matrix DX{pq) has eigenvalues with both positive and negative real 
parts. Assuming that k of these real parts are positive and the remaining n — k 
are negative, the stable manifold theorem ensures the existence of two invariant 
manifolds W'^{pq) and W^{po) with dimensions (\\m.W'^{po) = k and dimW^{po) = 
n — k, such that they intersect transversally one each other in po- 

Theorem 62 [1] Let pq be a nondegenerate strong singular point of the vector 
field X having an inverse Jacobi multiplier V defined in a neighborhood of po- 
Then V vanishes on W^{po) {resp. W^{po)) provided that divA'(po) < {resp. 
divA'(po) >0). 



X = 



\{-y + xf{x,y)), y = \{x + yf{x,y)) 



(35) 
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9.2 Time— dependent inverse integrating factors 

In [51], the authors consider autonomous second order differential equations 

X = w{x, x) , (36) 

with w G C°°{U) and W C an open set. They associated to (1361) the first order 
planar system defined on U in the usual way 



X 



y , y = w{x,y) . 



(37) 



Moreover, it is associated to equations fl36l) and fl37|) the vector fields X = dt + 
xdx + w{x, x)d± and X = yd^ + w(x, y)dy, respectively. A nonconstant function 
I{t,x,y) is called an invariant (or non-autonomous first integral) of system fl37|) 
in U if it is constant along the solutions of (1571) . In other words, XI = must be 
satisfied in U. Of course, we can find at most two functionally independent invari- 
ants of (jnZD- Notice that an invariant provides information about the asymptotic 
behavior of the orbits. 

A symmetry of (136!) is a diffeomorphism $ : (t, x) ^ (t, x) that maps the 
set of solutions of (!36|) into itself. Therefore, the symmetry condition for (!36|) 
is just x" = w{x,x'), where the prime denotes the derivative ' = d/dt. When 
the symmetry is a 1-parameter Lie group of point transformations then i = 



t + e^{t,x) + 0{e^ 



X 



X + er]{t, x) + 0(e ), for e close to zero, and the vector field 



y = ^(t, x)dt + r](t, x)dx is called the infinitesimal generator oi the 1-parameter Lie 
group of point transformations It is well known that the determining equations 
for Lie point symmetries can be obtained from the linearized condition 



y^'^\x — w{x,x)) = when x 



wix, X 



{31 



where 3^^^^ = 3^ + r/ W (t, x, x)d± + +77'^! (t, x, x, x)dx is the so-called second prolonga- 
tion of the infinitesimal generator y and rf^^ {t, x, x) = Dtrj — xDt$,, r/'^' (t, x, x, x) = 
DtT]^^^ — xDti where Dt = dt + xd^ + xd± is the operator total derivative with 
respect to t. Of course, since (l36l) is autonomous, it always admits the generator 
y = dt oi a. Lie point symmetry. Let Cr denote the set of all infinitesimal gener- 
ators of 1-parameter Lie groups of point symmetries of the differential equation 
(15^ . It is known that is a finite dimensional real Lie algebra, where we denote 
r = dim£j.. Moreover, for autonomous second order differential equation we have 
re {1,2,3,8}. 

For any = ^i(t, x)dt+r]i(t, x)dx G Cr, easily one can check that the Lie bracket 
[A", 3^]^'] = fii{t,x,x)X where fii{t,x,x) = X^i and yf^ = + rif\t,x,x)d± is the 
first prolongation of 3^. If r > 2, we define the functions 



v,,{t,x,x) = dei{x,yf\yf} 



1 X w{x,x) 

^i{t,x) rii{t,x) r]f\t,x,x) 
^j{t,x) Vj{t,x) r]f\t,x,x) 



(39) 
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for i,j G {1, . . . , r} with 1 < i < j < r. The aim of the work [51] is to generahze 
the concept of inverse integrating factor V{x,y) of system fl37l) via the functions 
Vij{t,x,y) defined in fl39l) . In fact, in the autonomous particular case dVij/dt = 0, 
we get that Vij is just an inverse integrating factor of fl37|) . On the contrary, when 
dVij/dt ^ 0, in [51] it is proved that the zero-sets V~^{0) and V^J^(O) have similar 
properties. The next result provides the connection between inverse integrating 
factors of system fl57|) and the functions Vij(t,x,y). 

Proposition 63 Assume that system (137|) possesses an r -dimensional Lie point 
symmetry algebra with r > 2 and define the functions Vij{t,x,x) as in (1391) . 

(i) Vij satisfies the linear partial differential equation XVij = Vy div^Y, where 
X = dt + xdx + w{x, x)dx- 

(ii) // r > 3 then, the ratio of any two nonzero Vij is either a constant or an 
invariant of (1371) . 

(iii) IfVj = 0, then {rji —yii)/{j]j — y^j) is an invariant of system ( 1571) . 

The next theorem is about the invariant curves of X contained in V^~^(0) and 
give us an extension of Theorem 9 in [63] for a case with dVij/ dt ^ 0. We put special 
emphasis on periodic orbits of (1571) of any kind (isolated and, therefore, limit cycles 
or non-isolated and so belonging to a period annulus). Recall here that a limit cycle 
7 := {{x{t),y{t)) eU : <t <T} is hyperbolic if §^diY X{x{t),y{t))dt ^ 0. On 
the other hand, a curve f{x, y) = defined on U is invariant for A? if Xf = Kf 
for some function K{x, y) called cof actor. 

Theorem 64 ^ LetU G be an open set and assume that x = w{x, x) with w 
smooth in U admits an r-dimensional Lie point symmetry algebra Cr with r > 2. 
Consider the functions Vij{t, x, x) defined in (I5U]) for i,j G {1, . . . , r} with 1 < i < 
j < r. Suppose that 7 = {x(t),y{t)) dU is a T -periodic orbit of (!37|) . Then the 
next statements hold: 

(i) If Vij{t,x,x) = V{x,x) ^ 0, with V G C^iU), then V{x,y) is an inverse 
integrating factor of system (1571) inU. In particular, if '-j is a limit cycle, 
then 7 C {V{x,y) = 0}. 

(ii) If Vij(t,x,x) = F{t)G{x, x) ^ with non-constants F and G G C^iU), then 
F = aF with a G M\{0} and G{x,y) = is an invariant curve of system 
(157|) . Moreover, we have: 

(ii.l) If ^ C {G = 0} and G is analytic on U, then G is not square-free, i.e., 
G{x,y) = g'^{x,y)u{x,y) with a positive integer n> 1 and g and u are 
analytic functions on lA satisfying 7 C {(7 = 0} and 7 ^ {u = 0}. 
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(ii.2) Ij ^ {G = 0} then 7 is hyperbolic and aT = f^divX{x{t),y{t))dt. 

An immediate consequence is obtained. 

Corollary 65 [54j Assume that x = w{x,x), with w smooth in the open set U C 
M^, admits an r -dimensional Lie point symmetry algebra Cr with r > 2. Consider 
the functions Vij(t,x,x) for i,j G {1, . . . ,r} with I < i < j < r. If there is one 
Vij{t,x,y) = F(t)G{x,y) ^ with non-constants F and G G C^iU), then system 
fl57]) does not have period annulus in U. 

In the sequel, we concentrate our attention in the 2-dimensionaI case £2- In [51] 
it is proved that, if 9f G £2, then the autonomous or separate time- variable forms 
of Vij{t, X, x) given in Theorem [61] are the only possibilities. Moreover, defining the 
domain of definition of the infinitesimal generators as the unbounded open strip 
E = {{t,x) G M X X} C M^, one has the following result. 

Theorem 66 |54j Assume that x = w{x,x) with w smooth in U <Z M? admits a 
2-dimensional Lie point symmetry algebra C2 spanned by the C^(S) vector fields 
yi = dt and 3^2 such that [3^1, 3^2] = Ci3^i + C23^2- 

(i) //c2 = and 3^2 ^ ('='); then Vi2{t,x,x) = G{x,x) with G{x,y) = y'^[ci + 
ya'{x) — 13' {x)] + (3{x)w{x,y) an inverse integrating factor of X in W = 
U n {X X M} provided that G ^0. Moreover, for analytic vector fields 3^2 iiT- 
H, X has no limit cycles in W . 

(ii) //c2 7^ then, changing the basis of C2 such that [3^1,3^2] = 3^i; we have that 
Vi2(t, x, x) = exp(c2t)G(x, x) with G{x, x) = x[c2xa{x) — C2l3{x) + x'^a'{x) — 
xl3'{x)] + P{x)w{x,x). In addition, dw/dx = or j3{x) = 0. If G ^ and 
U is a simply connected domain, then X has no periodic orbits in U and all 
the a or u) -limit sets of X are contained in the invariant curve G{x,y) = 
ofX. 

As an application of these results to polynomial Lienard systems, in [54j| it is 
proved the next theorem. 

Theorem 67 [54j The polynomial Lienard differential equation x + f{x)x + g{x) = 
with f,g& M[x] having a r-dimensional Lie point symmetry algebra Cr with r > 2 
has no limit cycles in M?. 

Acknowledgements 

We would like to thank Professor Hector Giacomini, from Universite de Tours 
(France), for his useful comments on this survey and for encouraging us to study 
the inverse integrating factor. 



52 



References 



[1] A. Andreev, Investigation on the behaviour of the integral curves of a system 
of two differential equations in the neighborhood of a singular point, Transla- 
tion Amer. Math. Soc. 8 (1958), 187-207. 

[2] Andronov, a. a. et al. Theory of Bifurcations of Dynamic Systems on 
a Plane, John Wiley and Sons, New York, 1973. 

[3] L.R. Berrone and H. Giacomini, On the vanishing set of inverse inte- 
grating factors, Qual. Th. Dyn. Systems 1 (2000), 211-230. 

[4] L.R. Berrone and H. Giacomini, Inverse Jacobi multipliers. Rend. Circ. 
Mat. Palermo (2) 52 (2003), 77-130. 

[5] T.R. Blows and L.M. Perko, Bifurcation of limit cycles from centers and 
separatrix cycles of planar analytic systems, SIAM Rev. 36 (1994) 341-376. 

[6] Bluman, G. W. and Anco, S. C. Symmetry and integration methods for 
differential equations. Applied Math. Sciences 154, 2002 (New York: Springer) 

[7] G. D. Bluman and S. Kumei, Symmetries and Differential Equations, 
Springer, New York, 1989. 

[8] L. Cairo, H. Giacomini and J. Llibre, Liouvillian first integrals for the 
planar Lotka-Volterra system. Rend. Circ. Mat. Palermo (2) 52 (2003), 389- 
418. 

[9] L. Cairo and J. Llibre, Integrability of the 2D Lotka-Volterra system via 
polynomial first integrals and polynomial inverse integrating factors, J. Phys. 
A 33 (2000), 2407-2417. 

[10] J. Chavarriga, Integrable systems in the plane with a center type linear 
part, Appl. Math. (Warsaw) 22 (1994), 285-309. 

[11] J. Chavarriga, A class of integrable polynomial vector fields. Appl. Math. 
(Warsaw) 23 (1995), 339-350. 

[12] J. Chavarriga and I.A. Garcia, The Poincare Problem in the Non- 
Resonant Case: An Algebraic Approach, Differ. Geom. Dyn. Syst. 8 (2006), 
54-68. 

[13] J. Chavarriga, I.A. Garcia and J. Cine, On the integrability of differen- 
tial equations defined by the sum of homogeneous vector fields with degenerate 
infinity, Internat. J. Bifur. Chaos Appl. Sci. Engrg. 11 (2001), 711-722. 



53 



[14] J. Chavarriga, B. Garcia, J. Llibre, J.S. Perez del Rio, J. A. 
Rodriguez, Polynomial first integrals of quadratic vector fields, J. Differ- 
ential Equations 230 (2006), 393-421. 

[15] J. Chavarriga, LA. Garcia and J. Sorolla, Resolution of the Poincare 
problem and nonexistence of algebraic limit cycles in family (I) of Chinese 
classification, Chaos Solitons Fractals 24 (2005), 491-499. 

[16] J. Chavarriga, LA. Garcia and J. Sorolla, Non-nested configuration 
of algebraic limit cycles in quadratic systems, J. Differential Equations 225 
(2006), 513-527. 

[17] J. Chavarriga, H. Giacomini and J. Cine, The null divergence factor, 
Publicacions Matematiques 41 (1997), 41-56. 

[18] J. Chavarriga, H. Giacomini and J. Cine, An improvement to Darboux 
integrability theorem for systems having a center. Appl. Math. Lett. 12 (1999), 
85-89. 

[19] J. Chavarriga, H. Giacomini and J. Cine, On a new type of bifurcation 
of limit cycles for planar cubic systems. Nonlinear Anal. 36, (1999), 139-149. 

[20] J. Chavarriga, H. Giacomini and J. Cine, Polynomial inverse integrat- 
ing factors, Ann. Differential Equations 16 (2000), 320-329. 

[21] J. Chavarriga, H. Giacomini, J. Cine and J. Llibre, On the integra- 
bility of two-dimensional flows, J. Differential Equations 157 (1999), 163-182. 

[22] J. Chavarriga, H. Giacomini, J. Cine and J. Llibre, Darboux integra- 
bility and the inverse integrating factor. J. Differential Equations 194 (2003), 
116-139. 

[23] J. Chavarriga, H. Giacomini and M. Grau, Necessary conditions for 
the existence of invariant algebraic curves for planar polynomial systems. Bull. 
Sci. Math., 129 (2005), 99-126. 

[24] J. Chavarriga and J. Cine, Integrability of a linear center perturbed by a 
fourth degree homogeneous polynomial. Publ. Mat. 40 (1996), 21-39. 

[25] J. Chavarriga and J. Cine, Integrability of a linear center perturbed by a 
fifth degree homogeneous polynomial. Publ. Mat. 41 (1997), 335-356. 

[26] J. Chavarriga and J. Cine, Integrable systems via inverse integrating 
factor, Extracta Math. 13 (1998), 41-60. 



54 



[27] J. Chavarriga, J. GiNE AND M. Grau, Integrable systems via polynomial 
inverse integrating factors, Bull. Sci. Math. 126 (2002), 315-331. 

[28] J. Chavarriga and M. Grau, A family of non-Darboux-integrable 
quadratic polynomial differential systems with algebraic solutions of arbitrarily 
high degree. Appl. Math. Lett. 16 (2003), 833-837. 

[29] J. Chavarriga, J. Llibre and J. Sotomayor, Algebraic solutions for 
polynomial systems with emphasis in the quadratic case. Exposition. Math. 15 
(1997), 161-173. 

[30] J. Chavarriga and M. Sabatini, A survey of isochronous centers. Qual. 
Theory Dyn. Syst. 1 (1999), 1-70. 

[31] C. Christopher, Invariant algebraic curves and conditions for a centre. 
Proc. Roy. Soc. Edinburgh Sect. A 124 (1994), 1209-1229. 

[32] C. Christopher, Liouvillian first integrals of second order polynomial dif- 
ferential systems., Electron. J. Differential Equations, Vol. 1999(1999), 1-7. 

[33] C. Christopher, J. Llibre, C. Pantazi and S. Walcher, Darboux 
integrating factors: Inverse problems. Preprint, 2008. 

[34] C. Christopher, J. Llibre and J.V. Pereira, Multiplicity of invariant 
algebraic curves in polynomial vector fields. Pacific J. Math. 229 (2007), 63- 
117. 

[35] A. Cohen, An Introduction to the Lie Theory of One-Parameter Groups, 
with Applications to the Solutions of Differential Equations, D. C. Heath, New 
York, 1911. 

[36] B. Coll, A. Ferragut and J. Llibre, Polynomial inverse integrating 
factors of quadratic differential systems. Preprint, 2008. 

[37] C.B. Collins, Conditions for a center in a simple class of cubic systems. 
Differential Integral Equations 10 (2) (1997), 333-356. 

[38] G. Darboux, Memoire sur les equations differentielles algebriques du premier 
ordre et du premier degre (Melanges), Bull. Sci. Math. 32 (1878), 60-96; 123- 
144; 151-200. 

[39] M.V. DoLOV, The structure of a single-valued integrating factor neaar a 
cycle. (Russian) Differentsial'nye Uravneniya 17 (1981), 1490-1492. 



55 



M.V. DOLOV AND S.A. Chistyakova, Algebraic differential equations with 
an integrating factor of Darboux type. (Russian) Izv. Akad. Nauk Respub. 
Moldova Mat. (1993), 96-106. 

M.V. DoLOV AND A. A. Alekseev, On the absence of limit cycles in dy- 
namical systems with an integrating factor of a special type, Differentsial'nye 
Uravneniya 30 (1994), 947-954; translation in Differential Equations 30 
(1994), 876-883. 

N. Du AND S.H. Chen, Formal integrating factors as a method of distin- 
guishing between a center and a focus. J. Math. (Wuhan) 17 (1997), 231-239. 

F. DuMORTlER, Singularities of vector fields on the plane. J. Differential 
Equations 23 (1977), 53-106. 

F. DUMORTIER, J. Llibre AND J.C. Artes, Qualitative theory of planar 
differential systems. Universitext, Springer- Verlag, Berlin Heidelberg, 2006. 

F. DUMORTIER, D. Panazzolo and R. Roussarie, More limit cycles than 
expected in Lienard equations. Proc. Amer. Math. Soc. 135 (2007), 1895-1904 
(electronic). 

F. DUMORTIER, R. Roussarie and C. Rousseau, Hubert's 16th problem 
for quadratic vector fields. J. Differential Equations 110 (1994), 86-133. 

A. Enciso and D. Peralta-Salas, Existence and vanishing set of inverse 
integrating factors for analytic vector fields, preprint, 2008. 

A. Ferragut and J. Llibre, On the remarkable values of the rational 
first integrals of polynomial vector fields. J. Differential Equations 241 (2007), 
399-417. 

A. Ferragut, J. Llibre and A. Mahdi, Polynomial inverse integrating 
factors for polynomial vector fields. Discrete Contin. Dynam. Systems, 17 
(2007) 387-395. 

LA. Garcia, H. Giacomini and M. Grau, The inverse integrat- 
ing factor and the Poincare map. Trans. Amer. Math. Soc, to appear. 
arXiv:0710.3238vl [math.DS] 

LA. Garcia, H. Giacomini and M. Grau, Generalized Hopf bifurca- 
tion for planar vector fields via the inverse integrating factor. Preprint, 2009. 
arXiv:0902.0681vl [math.DS] 

LA. Garcia and J. Gine, Generalized cofactors and nonlinear superposition 
principles, Appl. Math. Lett. 16 (2003), 1137-1141. 



56 



[53] LA. Garcia and J. Gine, Non-algebraic invariant curves for polynomial 
planar vector fields. Discrete Contin. Dyn. Syst. 10 (2004), 755-768. 

[54] LA. Garcia, J. Gine and S. Maza, Periodic Solutions of 2nd Order Dif- 
ferential Equations with 2-Dimensional Lie Point Symmetry Algebra, Preprint 
2009. 

[55] LA. Garcia, J. Gine and J. Sorolla, On the existence of polynomial 

inverse integrating factors in quadratic system,s with limit cycles, Dyn. Contin. 
Discrete Impuls. Syst. Ser. A Math. Anal. 12 (2005), 151-173. 

[56] LA. Garcia and D.S. Shafer, Integral invariants and limit sets of planar 
vector fields, J. Differential Equations 217 (2005), 363-376. 

[57] A. Gasull and H. Giacomini, A new criterion for controlling the number 
of limit cycles of some generalized Lienard equations, J. Differential Equations 
185 (2002), 54-73. 

[58] A. Gasull and H. Giacomini, Upper bounds for the number of limit cycles 
of some planar polynomial differential systems, Preprint 2008. 

[59] A. Gasull, H. Giacomini and J. Llibre, New criteria for the existence 
and non-existence of limit cycles in Lienard differential systems. Dynamical 
Systems: An International Journal, 2008 pp. 1-15, 2008. 

[60] H. Giacomini, J. Gine and M. Grau, Integrability of planar polynomial 
differential systems through linear differential equations. Rocky Mountain J. 
Math., 36 (2006), 457-485. 

[61] H. Giacomini, J. Gine and M. Grau, Linearizable planar differential 
systems via the inverse integrating factor, J. Phys. A: Math. Theor. 41 (2008), 
135-205. 

[62] H. Giacomini, J. Gine and J. Llibre, The problem of distinguishing 
between a center and a focus for nilpotent and degenerate analytic systems, J. 
Differential Equations 227 (2006), 406-426. 

[63] H. Giacomini, J. Llibre and M. Viano, On the nonexistence, existence, 
and uniqueness of limit cycles, Nonlincarity 9 (1996), 501-516. 

[64] H. Giacomini, J. Llibre and M. Viano, On the shape of limit cycles that 
bifurcate from Hamiltonian centers. Nonlinear Anal. 41, (2000), 523-537. 

[65] H. Giacomini, J. Llibre and M. Viano, The shape of limit cycles that 
bifurcate from non-Hamiltonian centers. Nonlinear Anal. 43 (2001), no. 7, Ser. 
A: Theory Methods, 837-859. 



57 



[66] H. GlACOMlNl AND M. Ndiaye, New sufficient conditions for a center and 
global phase portraits for polynomial systems. Publ. Mat. 40 (1996), 351-372. 

[67] H. GlACOMlNl AND M. VlANO, Determination of limit cycles for two- 
dimensional dynamical systems, Phys. Rev. E 52 (1995), 222-228. 

[68] H. GlACOMlNl, M. VlANO AND J. Llibre, Semistable limit cycles that 
bifurcate from centers, Internat. J. Bifur. Chaos Appl. Sci. Engrg. 13 (2003), 
3489-3498. 

[69] J. GiNE, The nondegenerate center problem and the inverse integrating factor. 
Bull. Sci. Math. 130 (2006), 152-161. 

[70] J. GiNE, On the centers of planar analytic differential systems, Internat. J. 
Bifur. Chaos Appl. Sci. Engrg. 17 (2007), 3061-3070. 

[71] F. Gonzalez-Gascon, Limit cycles and symmetries of dynamical systems, 
Phys. Lett. A 76 (1980), 205-208. 

[72] J. GUCKENHEIMER AND P. HOLMES, Nonlinear oscillations, dynamical sys- 
tems, and bifurcations of vector fields, in: Applied Mathematical Sciences, vol. 
42, Springer- Verlag, New York, 1986. 

[73] J. Hale and H. Koqak, Dynamics and bifurcations. Texts in Applied Math- 
ematics, 3. Springer- Verlag, New York, 1991. 

[74] M. HiRSCH, Differential Topology, New York: Springer-Verlag, 1976. 

[75] E.L. Ince, Ordinary Differential Equations, Dover, New York, 1956. 

[76] J. P. JoUANOLOU, Equations de Pfaff algebriques. Lecture Notes in Mathe- 
matics, 708. Springer, Berlin, 1979. 

[77] R.E. KooiJ AND C.J. Christopher, Algebraic invariant curves and the 
integrability of polynomial systems, Appl. Math. Lett. 6 (1993) 51-53. 

[78] J. Llibre, Integrability of polynomial differential systems. Handbook of dif- 
ferential equations, 437-532, Elsevier /North- Holland, Amsterdam, 2004. 

[79] J. Llibre, J. Medrado and P. da Silva, Limit cycles for singular per- 
turbation problems via inverse integrating factor. Bol. Soc. Parana. Mat., to 
appear. 

[80] J. Llibre and C. Pantazi, Polynomial differential systems having a given 
Darbouxian first integral. Bull. Sci. Math. 128 (2004), 775-788. 



58 



J. Llibre and G. Rodriguez, Configurations of limit cycles and planar 
polynomial vector fields. J. Differential Equations 198 (2004), 374-380. 

J. Llibre and X. Zhang, Darboux theory of integrability in taking into 
account the multiplicity. J. Differential Equations 246 (2009), 541-551. 

V.A. LUNKEVIGH AND K.S. SlBlRSKll, Integrals of a general quadratic dif- 
ferential system in cases of a center. Diff. Equations 18 (1982), 563-568. 

J.F. Mattei AND R. Moussu, Holonomie et integrales premieres, Ann. 
Sci. Ecole Norm. Sup. 13 (1980) 469-523. 

J. Moulin Ollagnier, About a conjecture on quadratic vector fields. J. Pure 
Appl. Algebra 165 (2001), 227-234. 

J. Moulin Ollagnier, Liouvillian integration of the Lotka-Volterra system., 
Qual. Theory Dyn. Syst. 2 (2001), 307-358. 

R. Moussu, Sur ('existence d'integrales premieres pour un germe de forme 
de Pfaff, Ann. Inst. Fourier (Grenoble) 26 (1976) 171-220. 

J.M. Page, Note on singular solutions. Am. J. of Math. XVIII, (1896), 
95-97. 

J.M. Page, Ordinary Differential Equations with an Introduction to Liess 
Theory of the Group of One Parameter, MacMillan, London, 1897. 

C. Pantazi, Inverse problems of the Darboux theory of integrability for pla- 
nar polynomial differential systems. Doctoral thesis, Universitat Autonoma de 
Barcelona, July 2004. 

D. Peralta-Salas, Note on a paper of J. Llibre and G. Rodriguez concern- 
ing algebraic limit cycles. J. Differential Equations 217 (2005), 249-256. 

H. PoiNCARE, Sur ^integration algebrique des equations differentielles du 
premier ordre et du premier degre, Rend. Circ. Mat. Palermo 11 (1897), 193- 
239. 

H. PoiNCARE, Les methodes nouvelles de la mecanique celeste. III, Paris, 
1899. 

M.J. Prelle AND M.F. Singer, Elementary first integrals of differential 
equations.. Trans. Amer. Math. Soc. 279 (1983), 215-229. 



59 



[95] G. Reeb, Sur certaines propietes topologiques des varietes feuilletees, in: 
W.T. Wu, G. Reeb (Eds.), Sur les espaces fibres et les varietes feuilletees, 
Tome XI, in: Actualites Sci. Indust., vol. 1183, Hermann et Cie, Paris, 1952. 

[96] R. ROUSSARIE, A note on finite cyclicity property and Hilbert's 16th problem. 
Dynamical systems, Valparaiso 1986, 161-168, Lecture Notes in Math., 1331, 
Springer, Berlin, 1988. 

[97] R. ROUSSARIE, Bifurcation of planar vector fields and Hilbert's sixteenth prob- 
lem. Progress in Mathematics, 164. Birkhauser Verlag, Basel, 1998. 

[98] J. A. Sanders and F. Verhulst, Averaging methods in nonlinear dynam- 
ical systems, in: Applied Mathematical Sci., vol. 59, Springer-Verlag, New 
York, 1985. 

[99] D. SCHLOMIUK, Algebraic particular integrals, integrability and the problem 
of center., Trans. Amer. Math. Soc. 338 (1993), 799-841. 

[100] D. SCHLOMlUK, Algebraic and geometric aspects of the theory of polynomial 
vector fields. Bifurcations and Periodic Orbits of Vector Fields, D. Schlomiuk 
(ed.), (1993), 429-467. 

[101] M.F. Singer, Liouvillian first integrals of differential equations., Trans. 
Amer. Math. Soc. 333 (1992), 673-688. 

[102] W.H. Steeb, Nonlinear autonomous dynamic systems, limit cycles and 
oneparameter groups of transformations. Letters in Math. Phys. 2, (1977), 
171-174. 

[103] M. ViANO, J. Llibre and H. Giacomini, Arbitrary order bifurcations 

for perturbed Hamiltonian planar systems via the reciprocal of an integrating 
factor, Nonhnear Anal. 48 (2002), Ser. A: Theory Methods, 117-136. 

[104] S. Walcher, On the Poincare problem, J. Differential Equations 166 (2000), 
51-78. 

[105] S. Walcher, Plane polynomial vector fields with prescribed invariant curves, 
Proc. Roy. Soc. Edinburgh Sect. A 130 (2000), 633-649. 

[106] S. Walcher, Local integrating factors, J. Lie Theory 13 (2003), 279-289. 

[107] C. E. WULFMAN, Limit cycles as invariant functions of Lie groups, J. Phys. 
A12, (1979), L73-L75. 

[108] Ye Yian-Qian et al., Theory of limit cycles, Translations of Math. Mono- 
graphs 66, Amer. Math. Soc, Providence, 1986. 



60 



[109] H. ZOLADEK. On algebraic solutions of algebraic Pfaff equations, Stud. 
Math. 114 (1995) 117-126. 

Addresses and e-mails: 

Departament de Matcmatica. Universitat de Lleida. 
Avda. Jaume II, 69. 25001 Lleida, SPAIN. 

E-mails: garcia@matematica.udl . cat, mtgrau@matematica.udl . cat 



61 



